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Results  of  an  extensive  study  of  the  modified  betatron  accelerators  ara 
reported.  The  topics  are  beam  injection,  equilibrium  and  stability.  -Three-  1 
dimensional  particle  simulations  of  Injections  were  performed.  The  results  show 
good  agreement  with  theory.  Various  types  of  equilibria  were  studied  using  'ifU 
simulations.  Results  are  compared  to  the  calculations  of  Finn  and  Manheimer. 
Beam  stability  Is  Investigated  using  a  combination  of  analytic  and  numerical 
techniques.  Growth  rates  of  the  negative  mass,  resistive  wall,  beam  breakup, 

Lend  Image  displacement  Instabilities  are  calculated.  A  new  theory  of  the  : 
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negative  mss  Instability*  keeping  toroidal  corrections  to  Maxwell's  equations. 

Is  derived.  Good  agreement  between  theory  and  simulations  Is  obtained.  The 
negative  mss  Instability  Is  shown  to  be  the  dominant  Instability.  Simulations 
show  that  It  can  disrupt  the  beam  In  a  time  of  order  one  microsecond.  Stabili¬ 
zation  may  require  a  beam  of  large  minor  radius.  Work  on  the  other  Instabilities 
show  that  measures  can  be  taken  to  make  their  effect  negligible. 
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I.  Introduction 


The  possibility  of  using  betatrons  to  accelerate  high  current  electron 
beans  to  high  energy  has  generated  renewed  Interest  In  these  devices.  If 
successful,  they  will  provide  a  source  of  high  power  beams  which  Is  much 
more  compact  and  economical  than  linear  machines.  In  order  to  overcome  the 
space  charge  repulsion  which  limits  current  In  the  conventional  betatron,1 
a  toroidal  magnetic  field  Is  added.  This  combination,  shown  In  Fig.  1,  Is 
the  Modified  Betatron  Accelerator  (MBA).  Successful  operation  of  the  modi¬ 
fied  betatron  requires  (a)  the  ability  to  Inject  a  high  current  (“10  kA) 
of  low  energy  (-2  MeV)  electrons  onto  a  cyclotron  orbit  In  the  presence 
of  a  strong  (1-5  kG)  toroidal  magnetic  field;  (b)  a  stable  beam  equilibrium 
for  the  duration  of  the  acceleration  period,  typically  thousands  of  revolu¬ 
tions;  and  (c)  extraction  of  the  beam,  which  has  been  accelerated  to  high 
energy  (-50  MeV),  preferably  In  a  time  equal  to  one  rotation  period.  If 
Instabilities  are  present,  their  cumulative  effect  must  be  small  enough 
such  that  macroscopic  disruption  of  the  beam  1$  minimal. 

During  the  period  of  our  contract  with  the  Office  of  Naval  Research 
(ONR)  (1  Sept  81  -  31  Oct  83),  we  concentrated  primarily  on  requirement 
(b),  with  some  work  being  done  on  (a).  This  choice  was  made  after  consid¬ 
eration  of  what  was  of  most  relevance  to  the  overall  modified  betatron  pro¬ 
gram  and  consideration  of  where  we  could  make  the  greatest  contributions. 
The  existence  of  a  stable,  or  at  any  rate  confined,  state  for  the  beam  over 
many  circulation  periods  Is  the  minimal  requirement  for  the  MBA.  In  addi¬ 
tion,  this  subject  can  be  addressed  with  confidence  using  a  combination  of 
analytical  and  numerical  techniques  which  serve  as  a  check  on  each  other. 
For  the  Issue  of  Injection,  one  must  rely  on  three-dimensional  numerical 
simulations  alone.  A  number  of  such  simulations  were  performed  under  this 
contract,  but  their  expense  was  such  that  a  detailed  study  of  injection  was 
not  feasible  with  the  computing  resources  available.  Concerning  the  Issue 
ef  beam  extraction.  It  was  felt  that  It  would  be  premature  to  start  work 
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on  this  subject,  given  the  uncertainty  relating  to  requirements  (a)  and 
(b).  Also,  there  are  several  uses  for  the  modified  betatron  which  would 
not  necessarily  require  extraction  of  the  beam,  e.g. ,  microwave  and  free- 
electron  laser  applications. 


This  report  describes  In  detail  the  work  we  have  done.  In  Chapter  II 
we  give  the  results  obtained  from  one-turn  Injection  simulations  of  the 
scheme  In  which  the  diode  is  placed  Inside  the  torus.  In  Chapter  III  we 
describe  Investigations  Into  different  types  of  beam  equilibria  and  how 
these  equilibria  are  affected  by  the  addition  of  a  spread  in  energy  or 
Increased  transverse  emittance  to  the  beam.  In  Chapter  IV  the  results  of 
extensive  analytic  and  numerical  Investigations  into  beam  stability  are 
presented.  The  main  emphasis  here  Is  on  the  negative  mass  instability,  but 
we  have  also  done  calculations  on  resistive  wall  and  beam  breakup  instabil¬ 
ities.  The  latter  occurs  In  devices  In  which  the  beam  passes  through 
acceleration  gaps.  In  the  modified  betatron,  injection  and  extraction 
ports,  and  slots  to  enhance  magnetic  flux  diffusion  Into  the  torus  can  give 
rise  to  the  same  type  of  Instability.  Chapter  V  gives  our  main  conclus¬ 
ions. 


II.  Beam  Injection 


Injection  of  a  high  current  beam  Into  the  modified  betatron  is  compli¬ 
cated  by  the  presence  of  the  toroidal  magnetic  field.  The  simplest  way  to 
overcome  the  problem  Is  to  place  the  diode  Inside  the  torus,  and  inject 
along  the  toroidal  field  lines.2  An  alternative  scheme  Is  to  have  the 
diode  outside  the  torus  and  to  Inject  along  temporarily  open  field  lines.3 
We  have  looked  only  at  the  first  scheme. 

In  the  modified  betatron  the  equilibrium  radial  position  of  the  beam 

centroid,  rc,  is  determined  by  equating  the  centrifugal  force  to  the 

2 

opposing  force  from  the  self  and  applied  field,  i.e.,  yV0  /r  = 

E*e1f  ♦  V0  (B*xt  +  B*elf),  where  y  is  the  beam  energy  and  Ve  Is  the  tor¬ 
oidal  velocity.  (We  are  using  units  In  which  m  *  c  *  1.)  As  shown  by 
Kapetanakos  et  al.2  there  Is  a  convenient  cancellation  of  the  space-charge 
depression  of  the  Injection  energy  y,  by  the  toroidal  correction  to  the 
self  magnetic  field  B*B  .  As  a  result,  the  equilibrium  radius  Is  given 

simply  by  rc  »  y.,8/B*xt  ,  where  6  *  (1-1/y2)1/2.  This  result  was  con¬ 
firmed  by  a  three-dimensional  PIC  simulation  using  IVORY.4  In  this  study, 
a  pulse  was  Injected  Into  a  torus  of  major  radius  r0  ■  1  m.  The  coordin¬ 
ate  system  used  was  cylindrical,  r,  e,  z,  where  the  z-axls  coincided  with 
the  major  axis  of  the  torus,  and  e  was  the  toroidal  angle  (Fig.  1).  Two 
cases  were  studied,  one  In  which  Injection  took  place  at  r  *  r0,  z  *  0,  e 
■  0  (centered-injectlon  case)  and  one  Ion  In  which  injection  was  at  r  * 
r0  ♦  4  cm,  z  ■  0,  6  ■  0  (offset -Injection  case).  The  Injected  pulse  had 
a  radius  of  2  cm  and  had  a  trapezoidal  time  profile  of  current  I  and  energy 
y.  Over  the  first  5  ns,  the  current  rose  from  1*0  to  10  kA,  and  energy 
went  from  y  *  2  to  5.  This  was  followed  by  a  10  ns  "body"  In  which  I  *  10 
kA,  y  *  5  and  finally  by  a  5  ns  tall  which  was  the  mirror  image  of  the 
head.  These  values  for  y  are  the  space-charge  depressed  values.  The 
Injection  energy  yj  was  larger:  y-j  *  7  In  the  body  of  the  pulse.  The 
applied  Be  and  Bz  fields  were  modeled  by  the  required  curl -free  forms.5 
The  values  of  Bz  and  the  external  field  Index  n  were  chosen  to  ensure 
that  the  main  body  of  the  beam  followed  a  stable  orbit  which  did  not 


Intersect  the  wall.  The  results  of  a  centered-injection  case  in  which  the 
beam  was  Initially  cold  are  shown  in  Fig.  2.  In  this  simulation,  we  chose 
n  *  0.25,  B2  (r  «  r0)  *  114  kG  and  B@  (r  «  r0)  *  3  kG.  When  com¬ 
bined  with  the  value  of  -y.,  this  yields  a  predicted  equilibrium  radius 
of  100  cm.  The  body  of  the  pulse  in  fact  propagated  quiescently  through 
the  first  turn  (which  is  as  far  as  the  simulation  ran)  at  exactly  this  rad¬ 
ius.  Figure  2  shows  that  the  head  of  the  pulse  is  drifted  into  the  upper 
wall  at  2  •  8.9  cm.  This  is  due  to  the  fact  that  in  the  head  (and  tail)  we 
did  not  have  radial  force  balance  because  of  the  lower  value  of  y.  As  a 
result,  there  was  a  net  radial  force  Fr  on  the  particles  causing  an  Fr 
*  Be  drift  in  the  z-dlrection.  The  resulting  deposition  of  energy  on  the 
wall  may  create  plasma  leading  to  vacuum  degradation  and  possibly  an  ion 
resonance  Instability  of  the  beam.6  However,  serious  damage  to  the  wall  is 
not  expected,  since  the  point  where  the  beam  contacts  the  wall  moves  along 
with  v  •  c.  The  area  of  contact  Is  about  r0ae  ar  =  30x4.5  cm2. 

In  the  offset  Injection  case  the  centroid  of  the  beam  was  injected  at 
r  ■  104  cm  and  the  applied  vertical  field  was  112  G.  The  injection  energy 
was  the  same  as  before,  so  that  the  theoretically  predicted  equilib¬ 
rium  radius  Is  rc  ■  105  cm.  As  seen  in  Fig.  3,  the  body  of  the  pulse 
propagated  around  the  torus  at  about  this  radius.  We  conclude  that  a  large 
radial  offset  in  the  Injection  does  not  adversely  affect  beam  propagation 
or  beam  quality. 


Figure  2.  Phase  space  plots  of  a  20  ns  beam  pulse  being  injected  at  8  * 
r  s  r0  =  100  cm,  z  =  0  Into  a  betatron.  Ibis  plot  was  made  1 
ns  after  the  beginning  of  Injection.  While  the  plots  are 
rectangular,  the  simulations  was  done  in  r-e  cylindrical  geo¬ 
metry.  Lengths  are  in  centimeters;  V0  is  normalized  to  c. 


111.  Beam  Equilibrium 


During  a  hypothetical  1  millisecond  acceleration  time  in  a  1  m  radius 
betatron,  the  beam  circulates  about  48,000  times.  It  is  therefore  essen¬ 
tial  to  have  accessible  self-consistent  equilibria  for  the  beam.  In 
analytic  and  numerical  work,  such  equilibria  form  the  starting  point  for 
the  study  of  beam  stability.  The  most  detailed  and  exact  calculations  of 
beam  equilibria  to  date  have  been  carried  out  by  Finn  and  Manheimer,7  using 
a  cold-fluid  model  of  the  beam.  Their  results  show  that  self-consistent 
equilibria  do  In  fact  exist  In  the  parameter  regime  of  the  proposed  NRL 
experiment  (10  kA,  2  kG  toroidal  field,  y  %  3).  These  equilibria  are 
analogous  to  the  slow  mode  of  rotation  of  a  cylindrical  laminar  beam  in  an 
axial  guide  field.  In  order  to  test  the  sensitivity  of  their  equilibria  to 
the  cold-fluid  assumption,  we  used  Finn's  code  EQUIL38  to  initialize 
IVORY  for  Yj  *  14,  I  *  10  kA,  Be  «  1  kG.  In  IVORY,  numerical  noise 
due  to  discrete  particle  effects  causes  deviations  from  the  perfect  lamin- 
arity  assumed  in  EQUIL3.  The  drift  surfaces  obtained  from  EQUIL3  are  shown 
In  Fig.  4.  The  fact  that  these  surfaces  are  closed  in  the  region  where  the 
beam  Indicates  that  an  equilibrium  exists.  In  the  simulation,  we  propa¬ 
gated  the  beam  for  twenty  revolutions  around  the  torus  during  which  time 
the  beam  rotated  twice  about  Its  own  axis.  No  drifting  off  center  was 
observed,  and  the  beam  minor  radius  remained  constant.  We  conclude  that 
the  EQUIL3  equilibria  are  not  significantly  affected  by  the  addition  of 
small  amounts  of  temperature. 

The  methods  of  Finn  and  Manheimer  break  down  as  the  transition  from 
diamagnetic  to  paramagnetic  rotation  of  the  beam  about  Its  own  axis  is 
approached.  The  transition  occurs  at  the  energy  at  which  the  self  forces 
of  the  beam  are  equal  to  the  weak  focusing  force  due  to  the  external  field 
Index.  For  a  uniform  density  beam  with  Budker's  parameter  v  and  beam  rad¬ 
ius  rt>»  the  transition  occurs  at 

>T  "  <4v  'oV>1/3 


(1) 


For  the  parameters  me  have  adhered  to  in  all  our  simulations,  namely  r^  = 

2  cm,  r0  *  1  m,  I  *  10  kA,  me  have  yj  ■  18.  Above  this  energy,  the 
toroidal  field  Is  not  needed  for  equilibrium.  However,  for  reasons  related 
to  beam  stability,  orbital  resonances  etc.,  the  toroidal  field  will  in 
practice  probably  be  maintained  throughout  the  acceleration  of  the  beam. 

No  detailed  cold-fluid  calculations  of  beam  equilibria  along  the  lines  of 
Ref.  7  have  been  carried  out  for  y  >  Yy.  In  our  simulations  we  have 
found  that  by  using  the  same  numerical  presciption  as  for  y  <  yt  me  can 
find  satisfactory  equilibria  above  transition. 

Since  betatron  equilibria  analogous  to  the  slow  rotation  mode  in  cold 
cylindrical  beams  have  been  found,  it  is  natural  to  ask  whether  other 
cylindrical  beam  equilibria  carry  over  to  the  betatron.  We  have  looked  at 
both  a  "fast  mode"  and  a  Kapchlnsky-Vladlmlrsky  (K-V)  distribution,  but 
neither  appears  to  a  modified  betatron  equilibrium,  at  least  in  the  regime 
me  examined.  In  the  "fast  mode"  the  beam  Is  again  laminar,  but  rotates  at 
approximately  the  cyclotron  frequency  oe/y  »  eBe/mcy,  mhere  Be  is  the 
toroidal  magnetic  field.  When  Initialized  In  this  manner  at  Be  *  1  kG,  y 
■  12,  the  beam  slowly  expanded  In  minor  radius,  and  1/3  of  the  particles 
mere  lost  to  the  mall  after  14  turns. 

If  a  beam  from  a  non-lmmersed  cathode  is  propagated  into  a  region  with 
an  axial  magnetic  field,  then  a  K-V  distribution  can  result.  The  beam 
rotates  at  half  the  cyclotron  frequency,  fle/2y,  and  transverse  emlttance 
supplies  the  additional  radial  force  needed  for  equilibrium.  With  Be  »  1 
kG  and  y  *  12,  a  beam  with  a  K-V  distribution  expanded  and  lost  over  half 
of  Its  particles  In  15  turns.  At  y  •  20  the  behavior  was  similar.  When 
the  toroidal  field  mas  turned  off  at  y  ■  20  {>  yt).  homever,  the  K-V 
distribution  shomed  no  significant  Increase  in  minor  radius.  At  present  we 
do  not  have  physical  explanations  for  these  observations. 

An  Issue  of  particular  Importance  with  regard  to  beam  stability  Is  the 
amount  of  energy  spread  and  transverse  emlttance  the  equilibrium  can  toler¬ 
ate.9  Lacking  self-consistent  models  to  Initialize  our  simulations,  we 
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decided  simply  to  add  these  effects  to  the  slow-mode  equilibrium  and  see 
how  the  beam  evolved.  Taking  the  equilibrium  at  y  =  12,  Be  *  1  kG,  we 
added  a  small  random  component  Ay  to  the  energy  of  each  particle,  drawing 

— /Zy  ^ 

from  a  Gaussian  of  the  form  e  ,  where  the  y is  the  width. 

We  find  that  the  beam  develops  spiral  arms,  made  up  of  particles"  evaporat¬ 
ing"  from  the  beam,  as  shown  In  Fig.  S.  These  particles  can  end  up  on  the 
wall.  If  -rth  Is  large  enough,  the  whole  beam  eventually  disappears  in 
this  manner.  This  behavior  can  be  understood  using  results  from  EQUIL3. 

Fo  the  equilibrium  shown  in  Fig.  4,  one  can  compute  the  drift-surfaces  of 
test-particles  with  energies  different  from  the  equilibrium  value.  For  by 
•  ±  0.1  (less  than  IX  deviation  from  y  *  12)  there  are  closed  diamagnetic 
drift-surfaces  only  for  particles  less  than  1  cm  from  the  beam  center,  as 
shown  In  Fig.  6.  For  by  •  ±  0.2,  the  diamagnetic  region  is  still  smaller, 
as  shown  In  Fig.  7.  Particles  outside  this  region  are  on  paramagnetic  sur¬ 
faces,  which  are  distorted  near  the  beam  and  give  the  spiraling  effect. 
Particles  with  energy  greater  than  the  equilibrium  value  form  one  arm  and 
those  with  less  energy  form  the  other  (see  Fig.  7).  The  arms  do  not 
rotate.  As  current  Is  lost  from  the  beam  through  these  arms,  the  diamag¬ 
netic  region  shrinks.  (In  the  low  current  limit  all  equilibria  are  para¬ 
magnetic.)  Thus,  one  can  have  an  "instability"  In  which  all  of  the  beam 
eventually  migrates  to  paramagnetic  surfaces,  which  may  Intersect  the 
wall.  Choosing  Yth  ■  0.5,  «  1  kG,  we  find  that  2/3  of  the  part¬ 

icles  are  lost  to  the  wall  after  about  four  circulations.  The  remaining 
particles  rotate  in  the  paramagnetic  direction,  and  the  energy  spread  is 
greatly  reduced.  For  yth  »  0.1  (Fig.  5),  less  than  10X  of  the  part¬ 
icles  are  lost  through  the  "arms",  and  the  remaining  beam  is  cold.  The 
ability  of  the  beam  to  sustain  a  temperature  spread  Improves  as  one  gets 
further  below  the  transition  energy.  For  y  •  4,  for  example,  Ref.  7  esti¬ 
mates  that  a  low  density  tall  with  by/y  *50X  can  be  confined  in  the 
beam. 

The  effect  of  Increasing  the  transverse  emlttance  of  the  beam,  leaving 
everything  else  fixed.  Is  what  one  would  expect:  the  beam  Is  mismatched 
and  expands  and  contracts  with  an  average  radius  larger  than  the  initial 
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Figure  S.  Plots  of  Minor  cross-section  of  beam  with  It  energy  spread  after 
(a)  4  and  (b)  8  turns  around  torus.  Arrows  show  direction  of 
Motion. 
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2  cm.  Emlttance  was  Introduced  by  giving  the  particles  random  perpendic¬ 
ular  momenta  Pj^  drawn  from  the  distribution  e"^i  ^Pth  where 

Pth  the  width.  A  picture  of  a  beam  with  y  *  12,  Be  *  1  kG,  p^/m 
*  0.28  Is  shown  In  Fig.  8.  The  beam  expands  and  contracts  between  radii  of 

3  and  2  cm.  The  addition  of  emlttance  creates  an  energy  spread,  and  the 
formation  of  spiral  arms  Is  observed  during  the  evolution  of  the  beam. 

These  are  not  very  pronounced,  however,  since  the  energy  spread  is  small. 
For  a  case  with  pth/m  *  0.4,  the  only  difference  noted  was  that  the  maxi¬ 
mum  beam  radius  was  about  4  cm.  No  particle  loss  was  observed  In  these 
simulations.  One  would  expect  that  one  could  obtain  a  matched  state  in  the 
presence  of  Increased  emlttance  by  making  the  beam  rotate  faster  about  its 
own  axis.  Recent  work  by  Grossmann,  Finn,  and  Manhelmer10  confirms  this. 


fi.  ^4, 
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IV.  Beam  Stability 


The  circulating  electron  ring  in  the  modified  betatron  possesses  a 
number  of  modes  of  oscillation  which  can  be  driven  unstable  through  inter¬ 
actions  with  each  other,  with  a  resistive  wall,  with  gaps  in  the  wall  of 
the  device,  etc.  These  instabilities  have  received  a  good  deal  of  attent¬ 
ion  In  connection  with  synchrotron  accelerators  and  storage  rings  (low  cur¬ 
rent,  high  energy).11  It  appears,  however,  that  for  the  high-current 
regime  we  are  concerned  with,  some  of  these  instabilities  are  qualitatively 
different  from  what  is  found  at  low  current.  The  self-fields  and  toroidal 
magnetic  field  play  important  roles. 

The  main  emphasis  in  our  work  has  been  on  the  negative  mass  instabil¬ 
ity,  since  it  appears  to  have  the  largest  growth  rate.  In  addition,  it  is 
the  Instability  most  readily  amenable  to  analytic  and  numerical  treatment. 
Our  main  tool  in  this  study  has  been  the  3-0  PIC  simulation  code  IVORY.  In 
conjunction  with  running  simulations,  we  have  reanalyzed  the  theory  of  the 
negative  mass  instability.  The  numerical  and  analytic  approaches  provide  a 
check  on  each  other.  For  the  resistive  wall  and  beam  breakup  instabili¬ 
ties,  we  have  performed  only  analytic  and  semi -analytic  studies  owing  to 
the  difficulty  of  performing  numerical  PIC  simulations. 

Since  there  Is  a  considerable  amount  of  algebra  associated  with 
derivations  for  the  above  Instabilities,  we  have  relegated  the  derivations 
to  the  appendices.  Only  the  relevant  physical  results  will  be  described 

here. 

A.  Description  of  Negative  Mass  Instability 

The  conventional  negative  mass  Instability  causes  a  toroidal  beam  to 
bunch  In  the  toroidal  direction.12  Experimentally,  It  was  found  to  limit 
the  intensity  of  the  beam  by  creating  a  spread  In  energy,  which  caused  the 
beam  to  spread  out.12  In  our  simulations  of  the  high  current  betatron, 
however,  very  little  toroidal  bunching  is  observed.  The  mode  we  see  is  a 
primarily  transverse,  almost  rigid  kinking  of  the  beam,  as  in  Fig.  9. 


Figure  9.  HI attrition  of  transverse  kinking  effect  of  negative  mass  Insta 
blllty.  Shown  ere  plots  (r  vs.  z  and  6  vs.  r)  of  the  beam  at 
successive  times  during  the  development  of  the  %  ■  1  negative 
mass  Instability.  In  the  r-z  plot,  particles  at  all  e-posltlons 
are  plotted.  Dimensions  are  In  cm. 


Some  description  of  our  simulation  scheme  is  necessary  in  order  to  under¬ 
stand  this  figure.  For  stability  simulations,  IVORY  uses  a  spatial  grid  to 
resolve  the  minor  cross-section  of  the  beam.  The  toroidal  dependence  of 
the  fields  Is  treated  by  a  sum  of  Fourier  modes.  In  most  of  our  simula¬ 
tions  to  date,  just  three  mode  numbers  are  resolved,  namely  0,±£  where  £  is 
an  Integer.  Thus,  only  a  section  of  2u/£  radians  of  the  torus  needs  to  be 
modeled.  The  other  £-1  sections  are  identical.  The  particles  are  pushed 
in  a  three-dimensional,  nonlinear  manner.  For  the  sake  of  econorry  in  com¬ 
puting,  just  nine  groups  of  particles  are  used  to  resolve  the  wavelength  we 
are  simulating.  Initially,  these  groups  are  identical  and  equally  spaced 
(Fig.  9a),  so  that  only  equilibrium  fields  are  produced.  We  then  initiate 
the  instability  by  giving  the  particles  a  small -amplitude  (less  than  1%  of 
the  beam  radius)  sinusoidal  displacement  from  their  equilibrium  position. 
The  instability  then  commences  to  grow  at  a  well-defined  linear  growth 
rate,  with  a  real  frequency  approximately  equal  to  £flz/y  where  nz  is 
the  cyclotron  frequency  associated  with  the  vertical  field  Bz.  The  per¬ 
turbed  fields  generally  grow  about  four  decades  in  energy  during  the  linear 
growth  regime.  As  their  energy  becomes  comparable  to  the  equilibrium  field 
energy  the  effect  of  the  Instability  on  the  ring  becomes  visible,  as  seen 
In  Fig.  9b.  Since  the  phase  velocity  of  the  unstable  wave  is  almost  equal 
to  the  particle  circulation  velocity,  the  sideways  displacement  of  each 
particle  increases  monotonically  in  time.  Particles  moving  to  larger  rad¬ 
ius  take  longer  to  go  around  the  torus,  and  so  slip  back  relative  to  those 
at  the  equilibrium  radius.  Conversely,  particles  moving  to  smaller  radii 
advance  relative  to  those  at  the  equilibrium  radius.  As  a  result,  the  beam 
eventually  breaks  In  the  toroidal  direction,  as  shown  in  Fig.  9c.  Soon 
after  this,  the  beam  strikes  the  wall  and  about  2/3  of  the  current  is 
lost.  The  remainder  continues  to  circulate  for  some  time  as  a  diffuse 
cloud.  For  our  typical  parameters  (10  kA,  6  MeV)  and  choosing  £*1,  this 
whole  process  takes  about  20  circulation  periods,  or  0.4  usee,  to  occur. 


B.  Linear  Growth  Rates  of  Negative  Mass  Instability: 

Theory  and  Simulations 

The  simulations  we  performed  of  the  negative  mass  instability  showed 
substantially  different  growth  rates  from  those  predicted  by  existing 
theories.5’13  These  theories  neglected  the  toroidal  corrections  to 
Maxwell's  equations  without  rigorous  justification.  By  including  these 
effects  to  first  order  in  the  parameter  e  =  a/r0  =  (minor  radius/major 
radius),  we  obtain  good  agreement  with  simulations  over  a  wide  range  of 
parameters.  The  details  of  this  calculation  are  in  Appendix  A.  A  compar¬ 
ison  between  the  theory  and  IVORY  simulations  for  different  £-numbers  is 
shown  in  Fig.  10.  Here  the  beam  parameters  are  y  -  12,  I  =  10  kA,  Be  =  1 
kG.  The  agreement  is  rather  good  even  up  to  £  =  20,  where  the  long  wave¬ 
length  approximation  £/r0  «  1/a,  used  in  the  theory,  is  starting  to  break 
down.  The  growth  rate  scales  approximately  as  t1/2.  It  is  noticeable 
that  the  agreement  for  £  =  1  is  poorer  than  for  £  =  2  to  8.  We  suspect 
that  some  second  order  toroidal  terms  dropped  from  the  theory  may  be 
responsible.  A  derivation  keeping  these  terms  is  planned  (see  Appendix 
A). 


The  variation  of  growth  rate  r  with  y  for  £  =  1  is  shown  in  Fig.  11. 
This  figure  contains  results  for  both  Bq  =  1  and  0  kG.  At  B@  =  1  kG, 
the  theory  predicts  a  sharp  cutoff  at  about  y  *  21.  Due  to  numerical  noise 
in  IVORY,  it  is  not  possible  to  see  zero  growth  rate.  One  can  only  put  an 
upper  bound  on  r.  At  y  =  25,  r  was  reduced  by  a  factor  of  four  from  the 
value  at  y  *  20.  For  Be  =  0  no  equilibrium  is  possible  for  y  <  18  as  we 
saw  In  Chapter  III.  The  theory  predicts  that  for  y  >  12,  there  is  no  nega¬ 
tive  mass  instability.  A  simulation  for  y  =  20  tends  to  confirm  this:  we 
find  r/c  <  2xl0"5  cm-1. 


son  between  theory  and  simulations. 
Bo  =  0  are  shown. 


Effects  of  Transverse  Emlttance  and  Enerc 
Instability 


Spread  on  Negative  Mass 


The  growth  rates  computed  for  the  negative  mass  instability  In  the 
last  section  are  very  large  in  practical  terms.  It  has  been  suggested5*9 
that  Landau  damping  due  to  a  spread  In  toroidal  rotation  frequencies  can 
contribute  significantly  to  the  stabilization  of  these  modes.  Since  the 
beams  we  are  concerned  with  have  V@  =  c,  a  spread  in  rotation  frequencies 
can  come  only  from  a  spread  Ar  in  the  radii  of  the  particle  orbits.  If  the 
particles  were  tied  to  the  field-lines,  then  the  toroidal  frequency  spread 
An  would  be 


For  the  general  case.  In  which  particles  are  not  tied  to  the  field  lines, 
Eq.  (2)  gives  an  upper  bound  on  the  toroidal  frequency  spread  to  first 
order  in  Ar/r0.  "Betatron''  oscillations  of  the  particles  across  field 
lines  can  add  corrections  of  order  (Ar/r0)2  to  Eq.  (2).  For  the  Finn- 
Manheimer  equilibria,  we  have  computed  Afi/fi  by  direct  measurement  of  part¬ 
icle  rotation  frequencies  in  our  simulations.  One  would  expect  to  achieve 
the  maximum  value  of  a n/n  at  the  transition  energy  (y-j-  8  18),  since 
poloidal  rotation  ceases  and  the  particles  follow  the  toroidal  field 
lines.  Equation  (2)  then  predicts  Afl/fl  -  4*10-2,  since  the  beam  diameter 
is  4  cm  and  r0  8  100  cm.  At  y  8  20  and  y  *  12  we  measured  Afl/fl  *  3xl0-2 
and  A n/n  *  6xl0-3  respectively.  The  deviations  from  Eq.  (2)  are  due  to 
poloidal  rotation. 


To  estimate  the  stabilizing  effect  of  the  frequency  spread,  we  did 
some  simulations  In  which  the  spread  was  zero.  This  was  accomplished  by 
replacing  the  particle  beam  with  a  ring  made  up  of  rigid  disks.  We  find 
that  for  y  ■  12,  the  growth  rates  obtained  are  not  significantly  different 
from  those  of  the  particle  equilibrium,  as  shown  in  Fig.  12.  For  i  *  4, 
the  equations  of  Ref.  9  would  predict  a  10%  reduction  in  the  growth  rate 
due  to  the  frequency  spread  which  Is  not  Inconsistent  with  our  results. 
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Figure  12.  Comparison  of  growth  rates  for  a  rigid  beam  model  with  growth 

rates  for  self-consistent  particle  equilibrium,  as  a  function  of 
toroidal  mode-number.  Also  shown  are  reduced  growth  rates  for  a 
beam  with  enhanced  energy  spread  or  enhanced  transverse  emit¬ 
tance. 


A  key  question  Is  whether  the  beam  equilibria  we  are  simulating  can  be 
modified  to  obtain  substantial  reductions  In  growth  rate.  In  Chapter  III, 
we  looked  at  two  ways  of  Increasing  the  spread  in  frequencies;  viz.,  by 
Introducing  an  energy  spread  on  the  beam,  and  by  increasing  the  transverse 
emlttance.  The  result  of  having  a  1%  energy  spread  is  to  decrease  the 
growth  rate  by  30%,  as  shown  In  Fig.  12.  We  saw  In  Chapter  III  however, 
that  further  Increases  In  energy  spread  lead  to  serious  loss  of  confine¬ 
ment.  By  Increasing  the  transverse  emittance  from  30  mrad-cm  to  90 
mrad-cm,  the  growth  rate  dropped  by  somewhat  less  than  30%  (Fig.  12). 

Again,  larger  decreases  can  be  obtained  only  at  the  expense  of  beam 
quality. 

0.  Resistive  Wall  Instability 

The  resistive  wall  Instability  Is  driven  by  the  phase  lag  between  the 
electric  and  magnetic  response  of  the  wall  Image  to  perturbations  of  the 
beam.  In  this  section,  we  summarize  the  results  on  this  instability 
obtained  in  Appendices  B  and  C.  There,  we  find  that  any  slow  mode  on  the 
beam  Is  driven  unstable  by  a  resistive  wall.  There  are  two  transverse  slow 
modes  with  frequencies  given  by  u  •  Wz/y  -  Qez/y  for  the  cyclotron 
mode  and  u  «  Mz/y  -  Wg  for  the  "E*B"  or  "space-charge"  mode.  Here 
ug  Is  the  slow  transverse  bounce  frequency  of  the  beam. 

For  the  slow  cyclotron  mode,  the  growth  rate  peaks  at  u  -  0,  which 
occurs  for  a  *  n%/nz.  Since  fi2  and  y  Increase  during  beam  accelera¬ 
tion,  this  constraint  is  satisfied  only  for  a  limited  time  period  At  which 
corresponds  to  a  band  width 

A  df  At 

^  3t 


Setting  Aw  equal  to  the  growth  rate  r,  we  obtain  the  expression  for  the 
number  of  e-fol dings  as  the  beam  goes  through  the  resonance: 


m  *  TT^Tn  *y73T  (4) 

Equation  (3)  Is  maximized  late  In  the  acceleration  cycle,  when  y  is  large 
and  the  resonant  value  £  Is  small. 

After  passing  through  u  -  0,  the  instability  transitions  to  a  lower 
growth  rate  regime  which  persists  for  the  remainder  of  the  acceleration 
cycle.  In  this  regime  the  number  of  e-foldings  is  given  by 

(m)eff  *  J  rdt  (*) 

where  the  Integration  Is  from  the  time  when  u  £  0  to  the  end  of  the  accel¬ 
eration. 

For  the  slow  space-charge  mode  there  is  no  u>  =  0  resonance,  because 
0  <  <*»B/(fi2/Y)  <  1.  Thus,  one  has  only  the  regime  of  adiabatic  growth 
corresponding  to  Eq.  (5).  This  regime  covers  the  entire  acceleration 
cycle. 


Besides  the  transverse  slow  modes,  there  is  also  a  longitudinal  slow 
space-charge  mode  which  Is  resistive-wall  unstable.  There  is  no  resonance 
associated  with  this  mode. 

He  Illustrate  the  relative  Importance  of  the  various  forms  of  the 
resistive  wall  Instability  for  the  typical  modified  betatron  parameters 
shown  In  Table  1.  A  stainless  steel  wall  with  normalized  conductivity 
4»o  »  5.24*106  cm”1  is  assumed.  For  the  cyclotron  mode  with  w  *  0,  we  find 
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TABLE  1.  TYPICAL  MODIFIED  BETATRON  PARAMETERS  USED  IN  EVALUATING 
RESISTIVE  WALL  INSTABILITY  GROWTH  RATES 


Toroidal  Magnetic  Field  B.  2.5  kc 

Vertical  Magnetic  Field  (Initial)  Bz  115  g 

Toroid  Major  Radius  Rq  100  cm 

Toroid  Minor  Radius  R  10  Cm 

Beam  Radius  a  1  cm 

Beam  Current  v  0.59 

Beam  Energy  y  7-100 

Acceleration  Time  t  3*107  cm 
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a  growth  rate  r  »  4*10"5  cm"1.  Equation  (4)  is  maximized  for  y  *  70,  for 
which  i  *  1.  Total  growth  is  rat  ■  3.6  near  u  *  0.  The  contribution  from 
Eq.  (5)  is  much  larger,  however,  giving  (TAt)eff  *  51. 

The  growth  rate  of  the  transverse  space-charge  mode  is  comparable  to 
the  cyclotron  mode,  but  growth  occurs  over  the  whole  acceleration  cycle, 
giving  (m)eff  *  150. 

The  preceding  analysis  ignores  any  thermal  spread  of  beam  energy, 
which  would  tend  to  damp  the  instabilities.  In  Appendix  C,  we  show  that 
for  a  spread  of  a  few  percent  In  the  initial  beam  energy,  damping  should  be 
significant.  We  saw  in  Chapter  III  however,  that  attempting  to  stabilize 
the  Instability  In  this  way  can  seriously  disrupt  the  beam  if  it  is  near 
the  transition  energy. 

Successful  operation  of  the  modified  betatron  requires  cutting  the 
growth  (rat)eff  to  about  unity.  Employing  more  highly  conducting  cavity 
walls  (e.g.,  copper)  can  reduce  growth  by  a  factor  of  six.  Reducing  the 
acceleration  time  by  an  order  of  magnitude  would  then  effectively  eliminate 
these  Instabilities.  Alternatively,  a  spread  of  a  few  percent  in  the  elec¬ 
tron  energy  may  be  sufficient  to  damp  out  the  modes. 

E.  Beam  Breakup  and  Image  Displacement  Instabilities 


Our  work  on  the  beam  breakup  and  image  displacement  Instabilities  was 
motivated  mainly  by  the  proposed  racetrack  induction  accelerator. 14  This 
device,  sketched  In  Fig.  13,  Is  a  recirculating  accelerator  like  the  modi¬ 
fied  betatron.  However,  It  uses  accelerating  gaps  instead  of  a  betatron 
field  to  Increase  the  beam  energy.  In  a  linear  device  the  beam  breakup 
Instability  arises  from  a  resonant  coupling  between  beam  transverse  oscil¬ 
lations  and  m*l  (m  »  azimuthal  mode  number)  electromagnetic  cavity  modes 
localized  to  the  acceleration  gaps.  This  coupling  can  result  in  large 
lateral  displacements  of  the  beam.15  The  instability  is  not  due  to  the 
accelerating  function  of  the  gap.  Any  discontinuity  in  the  wall  can  have  a 
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similar  effect.  Thus,  beam  Injection  and  extraction  ports  In  the  racetrack 
or  modified  betatron  accelerators  could  excite  beam  breakup.  Recirculating 
devices  differ  from  linear  devices  in  that  (1)  the  beam  keeps  passing  the 
same  gaps  and  (2)  there  could  be  a  resonant  interaction  between  the  beam- 
breakup  and  negative  mass  Instabilities.  As  we  saw  in  Section  D,  there  are 
two  transverse  slow  modes  and  one  longitudinal  slow  mode  on  the  beam. 

These  have  negative  energy  and  so  can  couple  unstably  to  the  gaps.  Coup¬ 
ling  to  the  longitudinal  mode  occurs  only  due  to  toroidal  curvature.  To 
model  the  effect  of  the  gaps,  we  used  a  gap  response  function  defined  by 


w  +  iw  oi/Q  -  to 
O  0 


where  Z_l  is  the  transverse  impedance  of  the  gap,  Q  is  the  quality  factor 
and  to0  Is  the  gap  resonant  frequency.  Details  of  the  calculations  are 
given  In  Appendices  D  and  E. 

In  Appendix  D  toroidal  terms  are  dropped  so  that  the  beam  behaves  as 
If  It  were  In  a  straight  but  periodic  system.  In  this  way  we  can  study  the 
effect  of  passing  by  the  same  gaps  repeatedly.  The  parameters  we  studied 
are  given  in  Table  2.  The  acceleration  gap  parameters  are  those  of  a  lin¬ 
ear  Induction  accelerator  developed  by  the  National  Bureau  of  Standards. 

We  find  a  beam  breakup  growth  rate  r  *  1.3*10-3  cm”1,  yielding  30  e-fold- 
Ings  during  the  acceleration.  This  amplification  can  be  cut  by  reducing 
Q.  If  a  value  of  Q  *  6  can  be  attained,  then  only  5.5  e-foldings  occur. 

In  order  to  do  more  exact  calculations  taking  account  of  beam  acceler¬ 
ation  and  transients,  we  used  the  code  BALTIC.16  Table  3  summarizes  thir¬ 
teen  runs,  varying  N  (the  number  of  gaps),  Q,  to0,  and  y,  but  keeping 
NZi/Q  fixed.  Cases  1-4  show  the  effect  of  changing  the  gap  resonant  fre¬ 
quency  to0.  The  Interaction  Is  greatest  when  wq  is  an  integral  multiple 
of  the  beam  circulation  frequency  n2/-r,  e.g.,  when  uq  -  0.17757  *  13 
02/y.  Cases  5  and  6  show  the  effect  of  dropping  y.  Case  7  shows  that  by 
decreasing  Q  from  60  to  10,  the  growth  rate  drops  significantly. 
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TABLE  2.  NOMINAL  RACETRACK  INDUCTION  ACCELERATOR  PARAMETERS 


Path  Lengths 

L  =  460  cm 

Drlfttube  Radius 

R  *  7  cm 

Beam  Radius 

a  =  1  cm 

Gui  Field 

Bz  *  2  kg 

(wc  =  1.173  cm-1) 

Beam  Current 

I  »  1  kA 

(v  =  0.0588) 

Beam  Energy 

U  *  0.4  -  40  MeV 

(y  =  1.5  -  80) 

Number  of  Revolutions 

50 

Number  of  Gaps 

N  =  4 

Acceleration  per  Gap 

AU  =  0.2  MeV 

• 

O 

II 

< 

Gap  Resonant  Frequency 

880  MHz 

Mode  Quality  Factor 

Q  -  60 

Gap  Transverse  Impedance 

15  ohms 

(2j./0  «  0.5) 

Gap  Width 

i  =  5  cm 

TABLE  3.  SUMMARY  OF  BEAM  BREAKUP  INSTABILITY  CALCULATIONS  WITH  BALTIC 


0 

Wq 

(cm"1) 

1 

r 

(10-4  cm 

60 

0.18 

80 

8.3 

60 

0.1732 

80 

3.3 

60 

0.17757 

80 

8.9 

60 

0.1787 

80 

8.8 

60 

0.18 

60 

10.0 

60 

0.17757 

60 

12.3 

10 

0.18 

80 

3.4 

60 

0.18 

80 

9.6 

60 

0.18 

1.5 

12.0 

60 

0.18 

1.5-80 

9.8 

10 

0.18 

1.5-80 

2.6 

6 

0.18 

1.5-80 

1.2 

6 

0.17757 

1.5-80 

1.6 

Cases  8  and  9  treat  two  evenly  spaced  gaps,  showing  that  the  growth 
rates  are  about  the  same  as  for  a  single  gap  (NZj_/Q  is  fixed).  In  cuses 
10-13  there  are  4  gaps  spaced  30  cm  apart,  and  the  beam  is  accelerated. 

The  average  growth  rate  is  shown.  From  these  we  see  that  with  Q  =  6  growth 
of  the  beam  breakup  instability  is  negligible  (<  3  e-foldings  occur). 

The  image  displacement  instability  arises  due  to  the  interruption  of 
the  beam  image  current  at  a  discontinuity  in  the  drift-tube  wall.  It  is 
the  long  wavelength  limit  of  the  beam-breakup  instability.  For  our  param¬ 
eters,  we  compute  a  peak  growth  of  9*10-5  cm”1,  and  this  occurs  only  over  a 
narrow  range  of  parameters.  Its  effect  is  therefore  negligible.  No 
evidence  of  coherent  growth  of  this  instability  was  seen  in  the  BALTIC 
runs. 

In  Appendix  E,  the  effect  of  the  accelerating  gaps  is  put  into  the 
full  toroidal  dispersion  relation.  In  that  way,  the  interaction  between 
the  beam  breakup  and  negative  mass  instabilities  can  be  investigated.  The 
strongest  interaction  would  be  expected  when  the  gap  resonant  frequency 
matched  a  harmonic  of  the  beam  circulation  frequency.  For  r0  =  70  cm, 
this  occurs  for  £  *  13.  The  main  result  of  this  work  is  that  the  coupling 
between  the  two  instabilities  is  weak.  No  strong  hybridization  occurs. 
Figure  14  shows  that  the  negative  mass  instability  is  dominant  and  its 
growth  rate  is  insensitive  to  the  presence  of  the  gaps  (Q  =  6  is  assumed). 


V.  Conclusions 


Our  work  on  the  modified  betatron  accelerator  has  dealt  with  a  range 
of  topics  involving  beam  injection,  equilibrium  and  stability.  For  the 
problem  of  injection,  we  performed  the  first  detailed  particle  simulations 
of  firing  a  beam  into  a  torus.  These  showed  that  the  body  of  the  injected 
pulse  propagates  in  a  well-behaved  manner  around  the  torus  at  the  theo¬ 
retically  predicted  radius.  The  head  and  tail  of  the  pulse,  which  are 
unmatched,  drift  rapidly  in  the  vertical  direction,  striking  the  wall.  No 
serious  wall  damage  is  expected  to  result. 

On  the  subject  of  beam  equilibrium,  we  attempted  to  find  in  the  beta¬ 
tron  the  counterparts  of  beam  equilibria  found  in  linear  devices.  Follow¬ 
ing  the  work  of  Manheimer  and  Finn,  we  sucessfully  simulated  the  slow  rota¬ 
tion  mode  of  the  beam.  We  found  that  this  type  of  equilibrium  was  rela¬ 
tively  Insensitive  to  the  addition  of  small  amounts  of  transverse  temper¬ 
ature.  For  the  fast  mode  of  rotation  and  for  a  Kapchinsky-Vladimirsky 
distribution,  no  equilibrium  could  be  found  for  the  parameters  we  chose. 

We  emphasize  that  our  search  for  these  alternative  equilibria  was  not 
exhaustive. 

The  addition  of  an  energy  spread  or  of  increased  transverse  emittance 
to  the  beam  is  commonly  advocated  in  order  to  improve  beam  stability.  We 
Investigated  the  effect  these  changes  have  on  the  slow-mode  equilibrium. 

We  found  that  unless  the  beam  energy  is  far  removed  from  the  diamagnetic- 
paramagnetic  transition  energy,  the  tolerance  of  the  beam  for  an  energy 
spread  is  very  limited.  The  beam  tends  to  develop  spiral  arms,  along  which 
the  mismatched  particles  escape.  Significant  loss  of  particles  can 
result.  The  development  of  the  arms  is  explained  using  the  Finn-Manheimer 
model.  By  adding  transverse  emittance  to  the  beam,  large  oscillations  of 
the  minor  beam  radius  could  be  produced.  No  significant  loss  of  confine¬ 
ment  resulted. 


The  bulk  of  our  work  concerned  beam  stability,  and  the  negative  mass 
Instability  In  particular.  Analytic  and  numerical  methods  were  used.  A 
new  analytic  dispersion  relation  was  derived  keeping  toroidal  corrections 
to  Maxwell's  equations.  The  new  theory  overcomes  some  intrinsic  deficien¬ 
cies  of  previous  theories.  He  have  compared  growth  rates  for  the  negative 
mass  Instability  obtained  from  the  theory  with  3-0  numerical  simulations 
using  IVORY.  Good  agreement  is  obtained  over  a  wide  range  of  parameters. 
The  growth  rate  of  the  negative  mass  instability  is  so  large  that  it  can 
seriously  disrupt  the  beam  in  a  microsecond.  Addition  of  an  energy  spread 
and  transverse  emittance  give  reductions  on  the  order  of  30%  in  growth 
rate.  Much  larger  reductions  clearly  are  needed.  Based  on  our  work  to 
date,  we  believe  that  this  can  be  accomplished  only  by  starting  with  a  beam 
whose  radius  is  comparable  to  the  minor  radius  of  the  torus. 

In  addition  to  the  negative  mass  instability  we  looked  at  resistive 
wall,  beam  breakup  and  image  displacement  instabilities.  In  general,  these 
have  growth  rates  substantially  less  than  the  negative  mass  instability. 

The  resistive  wall  instability  can  probably  be  eliminated  by  using  a  copper 
lining  for  the  torus  and  decreasing  the  acceleration  time  to  hundreds  of 
microseconds.  The  beam  breakup  and  image  displacement  instabilities  were 
considered  for  their  relevance  to  recirculating  devices  which  use  accelera¬ 
tion  gaps.  He  found  that  by  using  low  qual ity-f actor  gaps  the  growth  can 
be  made  negligible. 
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I.  INTRODUCTION 


An  accurate  cold  beam  linear  dispersion  relation  for  the  negative  mass 
Instability  In  high  current  conventional  or  modified  betatrons  Is  critical 
for  designing  experiments,  benchmarking  computer  simulations,  and  develop¬ 
ing  Marm  beam  dispersion  relations.  In  this  report  we  derive  a  dispersion 
relation  correct  to  first  order  In  the  toroidal  aspect  ratio  of  the  beta¬ 
tron  valid  for  small  toroidal  mode  numbers.  Limited  comparisons  with  ear¬ 
lier  work  are  provided  for  reference. 

The  two  earlier  treatments  of  negative  mass  instability  growth  in  high 
current  betatrons  Included  toroidal  curvature  effects  In  the  particle 
dynamics  but  not  In  the  electromagnetic  field  equations.1*2  In  other 
words,  the  dispersion  relations  were  derived  for  beam  motion  in  a  toroidal 
cavity  but  with  fields  computed  for  an  off-center  beam  In  a  straight  tube 
of  circular  crossection.  The  first  analysis,  that  of  Sprangle  and 
Vomvoridls,  used  an  equation  of  continuity  appropriate  to  a  straight  cyl¬ 
inder,  as  well.  In  contrast,  the  second  analysis,  due  to  Hughes  and 
Godfrey,  Incorporated  an  equation  of  continuity  correct  for  a  toroidal 
beam.  Predicted  growth  rates  from  these  two  models  often  have  been  fc 
to  disagree  by  factors  of  two  for  parameters  corresponding  to  peak  growth 
and  by  even  more  at  higher  beam  energies.  More  disturbingly,  both  models 
systematically  overestimate  Instability  growth  rates  observed  in  three- 
dimensional  computer  simulations.3 »4  We  attribute  these  discrepancies 
to  Incomplete  treatment  of  toroidal  curvature  effects. 


The  dominant  curvature  correction  to  the  field  equations  is  easily 
Identified.  Radial  motion  of  an  electron  ring  (see  Fig.  1)  gives  rise  to 
azimuthal  bunching.  This  constitutes  the  difference  between  charge  contin¬ 
uity  In  cylindrical  and  toroidal  geometry.  Conversely,  azimuthal  bunching 
gives  rise  to  a  radial  electric  field,  which  constitutes  the  (main)  differ¬ 
ence  between  electromagnetic  fields  In  cylindrical  and  toroidal  geometry. 


1 


Figure  1.  Schematic  representation  of  betatron  cavity  with  electron  beam 
Mirror  and  optional  toroidal  magnetic  fields  are  applied. 
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The  additional  coupling  between  radial  and  azimuthal  beam  motion  need  not 
be  small*  Indeed,  for  sufficiently  large  energies  this  coupling  through 
the  fields  must  dominate  coupling  through  the  particle  dynamics.  The 
latter  falls  off  with  energy,  while  the  former  does  not. 

The  present  paper  Includes  these  and  other  (lesser)  curvature  effects 
to  first  order  In  the  equilibrium  and  perturbed  fields,  producing  what  we 
believe  to  be  a  more  accurate  dispersion  relation.  Like  its  predecessors, 
this  analysis  Ignores  the  Internal  cfynamics  of  the  beam,  treating  it  as  a 
string  of  rigid  disks.  The  assumption  of  long  azimuthal  wavelengths  also 
is  Implicit  In  our  derivation.  Relaxing  the  last  constraint  is  conceptu¬ 
ally  straightforward  but  algebraically  difficult. 

In  Section  II  we  assemble  the  various  particle  and  field  equations 
Into  a  5x5  self-adjoint  matrix  operator  acting  on  the  perturbed  beam  cen¬ 
troid  position  and  the  electrostatic  and  electromagnetic  potentials.  The 
differential  equations  are  solved  to  first  order  In  the  toroidal  aspect 
ratio  In  Section  III,  and  the  remaining  algebraic  equations  collapsed  to 
the  desired  dispersion  relation.  When  no  toroidal  magnetic  field  exists 
(l.e.,  for  a  conventional  betatron),  an  analytic  growth  rate  expression  for 
the  negative  mass  Instability,  valid  over  a  wide  parameter  range,  can  be 
derived  from  the  general  dispersion  relation,  and  this  is  done  in  Section 
IV.  We  have  not  yet  attempted  to  obtain  a  corresponding  simple  growth  rate 
expression  for  the  modified  betatron. 

Finally,  Section  V  presents  a  reevaluation  of  negative  mass  growth  for 
the  Office  of  Naval  Research  racetrack  induction  accelerator  design,5  based 
on  the  new  dispersion  relation.  The  average  instability  growth  rate  during 
acceleration  from  a  few  to  50  MeV  Is  reduced  by  a  factor  of  two  relative  to 
an  earlier  prediction.6  With  this  Improvement,  the  proposed  accelerator 
may  be  able  to  "outrun"  the  Instability  with  only  a  modest  beam  thermal 
energy  spread.7 

Comparisons  with  computer  simulation  results,  as  applied  to  the  Naval 
Research  Laboratory  modified  betatron  design,  are  described  elsewhere.8 
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II.  FORMULATION  OF  THE  PROBLEM 


We  consider  an  electron  ring  circulating  in  an  azimuthally  symmetric 
conducting  cavity.  Dimensions  and  coordinates  are  depicted  In  Figure  1. 

For  the  present  we  assume  only  that  the  beam  and  cavity  crossectlons  are 
symmetric  about  z«0.  Otherwise,  the  crossectlons  may  assume  any  reasonable 
shape.  The  characteristic  minor  radius  of  the  beam  must  be  much  less  than 
that  of  the  cavity  so  that  the  Internal  dynamics  of  the  beam  can  be 
Ignored.  A  mirror  magnetic  field  maintains  the  beam  In  a  circular  orbit, 
and  a  toroidal  field  can  be  added  to  confine  the  beam  against  Its  self¬ 
fields,  which  need  not  be  small. 

Linearized  equations  for  small  transverse  displacements  of  the  beam 
centroid  are  easily  obtained  from  the  single  particle  equations  of  motion. 


Perturbed  quantities  are  preceded  by  a  delta  (e.g.,  6z),  while  unperturbed 
quantities  are  not.  Total  time  derivatives  of  perturbed  particle  quanti¬ 
ties  are  represented  by  dots  above  the  quantities  (e.g.  6z).  Note  that  the 
last  term  In  Eq.  (2)  can  be  rearranged.  If  desired,  by  means  of  the 
equilibrium  radial  force  balance  equation, 

EP  ♦  v,  4  i  ve2'R  ’  0  <3> 


Although  the  perturbed  azimuthal  velocity  equation  can  be  derived  in 
many  ways,  using  the  single  particle  energy  formula  seems  simplest. 

E. 
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The  perturbed  azimuthal  angle  66  of  a  beam  disk  Is  related  to  6Ve  by 

ve 


6V0  «  R  66  +  -£-  6r 


(5) 


Perturbed  beam  currents  resulting  from  the  centroid  displacements 
described  by  Eq.  (1),  (2),  (4),  and  (5)  are 


6J. 


p  6r,  6J 


p  6z 


and 
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(6) 


Perturbed  charge  density  Is,  In  turn,  derived  by  substituting  the  perturbed 
currents  Into  the  continuity  equation  and  Integrating  the  result  in  time. 


fip  +  P«e  ♦  7  If  r  P6r  +  P6z  «  0 


(7) 


Alternatively,  Eq.  (7)  can  be  obtained  by  considering  how  the  density  of  a 
beam  element  changes  as  It  Is  displaced  Infinitesimally  in  each  direction. 

As  explained  In  the  Introduction,  we  are  limiting  consideration  to  low 
frequencies  and  long  azimuthal  wavelengths.  In  this  limit  we  need  only 
determine  the  electrostatic  and  azimuthal  electromagnetic  potentials,  which 

satisfy 
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The  electric  and  magnetic  field  combinations  appearing  In  Eq.  (1),  (2),  and 
(4)  are  expressed  In  terms  of  these  potentials  as 
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These  equations  complete  our  model. 

Let  us  now  cast  the  equations  In  matrix  form  for  compactness  and  to 
emphasize  their  symmetry.  In  so  doing  we  also  Fourier  transform  the  equa¬ 
tions  in  time  and  azimuthal  angle,  i.e.. 
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and  eliminate  6Ve  and  the  perturbed  electric  and  magnetic  fields. 
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Note  that  the  matrix  operator  is  self-adjoint,  as  one  would  hope.  Certain 
symbols  appearing  in  Eq.  (13)  are  defined  as 
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Equation  (13)  can  be  solved  numerically  for  beams  and  cavities  of 
arbitrary  cross  section.  (Before  doing  so.  It  probably  would  be  desirable 
to  perform  a  conformal  transformation  on  r,2  to  map  the  beam  and  cavity 
boundaries  onto  coordinate  surfaces.)  In  the  next  Section  we  Instead 
assume  the  beam  and  cavity  to  be  circular  and  concentric  In  cross  section, 
as  Illustrated  In  Fig.  1,  In  order  to  obtain  an  analytical  dispersion 
relation.  We  should  bear  In  mind,  however,  that  other  configurations  may 
exhibit  Improved  stability. 


III.  DERIVATION  OF  THE  DISPERSION  RELATION 


To  develop  a  dispersion  relation  from  Eq.  (13),  even  for  the  simple 
geometry  explicitly  represented  by  Fig.  1,  is  a  long  and  involved  process. 
For  the  sake  of  brevity,  we  here  restrict  ourselves  to  outlining  the  pro¬ 
cedure  and  citing  pertinent  results. 

Converting  the  partial  differential  operators  in  Eq.  (13)  to  a  set  of 
coupled  ordinary  differential  operators  is  the  first  step.  To  do  this  we 
carry  out  the  coordinate  transformation 

z  *  x  cos<j>,  r  *  R  +  x  sinY 

and  expand  the  potentials  as  Fourier  series  in  poloidal  angle  y. 

<«  -  £«♦"<*>  e1m*.  6A  »  i6Am(x) 

Next,  we  recast  Eq.  (13)  as  a  variational  integral,  insert  the  6$  and  6A 
expansions,  explicitly  Integrate  over  Y.  and  perform  variations  with 
respect  to  6Ym  and  6Am.  This  results  in  an  infinite  set  of  different¬ 
ial  equations  in  x^  coupling  all  the  poloidal  modes  of  the  potentials. 

Expanding  the  differential  equations  in  x  with  respect  to  the  cavity 
aspect  ratio,  b/R,  leads  to  a  natural  truncation  of  the  infinite  system. 
Basically,  an  expansion  to  order  |m|  in  b/R  is  consistent  with  dropping  all 
higher  poloidal  mode  numbers.  Two  orderings  spring  to  mind.  In  the  first 
<»,  t/ R  and  the  electron  oscillation  frequencies  all  are  taken  to  be  of 
order  unity  (or  smaller),  and  the  equations  expanded  to  order  b/R.  In  the 
second,  w  and  l/R  are  taken  to  be  of  order  b/R,  and  the  equations  expanded 
to  b2/R2.  The  first  alternative  Is  much  simpler,  seems  to  capture  the 
essential  physics,  and  agrees  reasonably  well  with  simulation  results,  so 
we  use  It  In  the  following  calculations.  We  are,  however.  Investigating 
the  consequences  of  the  second  ordering  and  will  publish  our  findings  at  a 
later  date. 


With  the  ordering  selected,  straightforward  but  tedious  calculations 
yield  for  x  <  a: 
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The  validity  of  these  solutions  to  the  truncated  equations  has  been  veri¬ 
fied  using  the  symbolic  manipulation  program  MAXIMA. 

Back  substitution  of  the  potentials  Into  the  variational  Integral 
allows  the  remaining  Integrations  to  be  carried  out  explicitly.  (For  con¬ 
sistency,  the  results  are  truncated  to  first  order  In  b/R.)  Again,  MAXIMA 
was  used  to  check  the  manipulations.  The  resulting  matrix  equations  for 
the  perturbed  centroid  location  appear  as  Eq.  (20).  Its  determinant  is 
the  sought-for  dispersion  relation. 
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In  evaluating  the  matrix  elements.  It  Is  necessary  to  know  the  equi¬ 
librium  fields.  They  are9 

Er  *  7(r-R)+W  (^7  +4£na)  (21) 

2  .  2 

B2  *  -»,  f  tr-R)  .  V,  *L-  (-  iy  ♦♦♦**»!) 

♦  >2°  -  "  T5!  (22) 

with  similar  expressions  for  Ez  and  Br.  Eq.  (3)  determines  the  magni¬ 
tude  of  B2°,  and  Be°  Is  arbitrary. 

The  dispersion  relation  Is 

(f  '  “z2)  (f  *  “r2  -  x/e)  -  u?  Be2/^  "  0  (23) 

It  represents  two  longitudinal  (m  ■  0)  modes,  described  by  the  longitudinal 
dielectric  constant, 

*  s  °2  -  7  (i  * *  *"  ?)  (7  -  “2)  (24) 

and  four  transverse  (m  *  ±1)  modes,  described  by  Eq.  (33)  with  x  *  0.  The 
transverse  oscillation  frequencies  uz  and  up  are. 


“r  5  "  (1“">  veT'^  ’  t*  ’  W  (26) 

Er  and  Bz  In  Eq.  (26)  are  given  by  Eq.  (21)  and  (22)  with  r«R  and 
Bz°  r-  itted.  The  last  term  In  Eq.  (26)  Is  very  small.  Typically,  the 
external  field  Index  n  Is  chosen  such  that  wr2  »  u^2.  Note  that  v  is 
Budker's  parameter,  equal  to  pa2/4. 

The  key  result  of  our  analysis  Is  the  coupling  coefficient  between 
longitudinal  and  transverse  modes,  which  determines  the  negative  mass 
Instability  growth  rate. 


For  comparison,  the  dispersion  relations  from  Ref.  1  and  2  also  can  be  cast 
In  the  form  of  Eq.  (23)  but  with  coupling  coefficients,  respectively. 


It  should  be  noted  that  there  Is  a  degree  of  arbitrariness  In  the 
choice  for  the  functional  form  of  p(r,z).  In  the  preceding  analysis  we 
assumed  p  to  be  constant  out  to  the  beam  minor  radius  a,  where  It  drops 
abruptly  to  zero.  One  might  Instead  have  chosen  pr/R  to  be  constant,  for 
instance.9  Fortunately,  such  changes  lead  only  to  Insignificant  modifica¬ 
tions  of  x.  and  u  2. 


IV.  Be°  -  o  GROWTH  RATE  FORMULA 


In  the  absence  of  a  toroidal  magnetic  field,  the  negative  mass  Insta¬ 
bility  dispersion  relation  reduces  to 

n2  -  «r2  -  x/£  -  0  (30) 

A  simple  analytical  growth  rate  expression  can  be  derived  from  Eq.  (30)  In 
a  reasonably  straightforward  manner,  If  n2  Is  much  less  than  wr2  and  can 
be  dropped  from  the  equation.  We  address  the  validity  of  this  assumption 
at  the  end  of  this  Section. 

Solely  for  the  sake  of  algebraic  simplicity,  we  invoke  four  approxima 

tlons. 


All  four  are  well  satisfied  for  typical  cases  of  Interest.  Then, 
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Significant  cancellations  occur  when  Eq.  (32)  Is  expanded, 
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The  term  linear  In  n,  In  fact,  only  survives  due  to  the  difference  In  the 


factors 
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which  appear  In  the  first  and  second  terms,  respectively,  of  the  definition 
for  x.  Modifications  to  the  expression  for  x  discussed  near  the  end  of  the 
preceding  Section,  due  to  changing  the  functional  form  of  p,  are  of  higher 
order  than  the  terms  In  Eq.  (33)  and  so  drop  out.  Nonetheless,  we  should 
not  be  surprised  If  more  accurate  dispersion  relation  derivations,  perhaps 
Including  high  frequency  electromagnetic  effects,  adjust  the  linear  term 
somewhat. 


Substitution  of  Eq.  (31)  and  (33)  Into 
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leads  Immediately  to  the  desired  growth  rate  formula. 
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In  obtaining  Eq.  (34),  we  have  taken 
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which  Is  accurate  for  y2  not  too  small,  an  assumption  we  have  already 
made.  Equation  (34)  reproduces  the  exact  numerical  solutions  of  the  dis¬ 
persion  relation  presented  In  the  next  Section  to  an  accuracy  of  10%. 

Equation  (34)  predicts  an  Instability  cutoff  for 

vy  >  4(1  -  n)  ^  +  2m  (35) 

At  smaller  values  of  vy,  the  negative  mass  instability  reduces  to  the  well 
known  expression10*11 

r  /  \*11/2 

r  "  i  [Tn  ^5  (l  +  2in  l)j  (36) 

Me  remark  for  completeness  that  Eq.  (28),  from  Ref.  1,  and  Eq.  (29),  from 
Ref.  2,  also  yield  high  vy  cutoffs.  In  the  case  of  Eq.  (29)  Instability 
ceases  for2*12 

vy  >  4(1  -  n) / ( 1  +  2in  j)  (37) 

while  Eq.  (28)  has  the  same  limit  without  the  factor  of  four.  Both  analy¬ 
ses  lead  to  Eq.  (36)  for  sufficiently  small  vy. 

Let  us  return  to  our  Initial  assumption. 


Based  on  Eq.  (34),  this  Inequality  becomes 

(nr)  7  (? +  2ln  l)  €  1 


or  for  typical  parameters 
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He  see  that  Eq.  (34)  Is  quite  generally  accurate  for  Be°=  0. 
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V.  THE  ONR  RACETRACK  INDUCTION  ACCELERATOR 


The  Office  of  Naval  Research  together  with  the  Naval  Research  Labora¬ 
tory  have  developed  the  design  for  a  racetrack  induction  accelerator 
Intended  to  accelerate  a  1  kA  (v  «  0.0588)  electron  beam  to  40  MeV  {y  =  80) 
in  fifty  revolutions.  In  an  earlier  report,6  we  Investigated  the  beam 
breakup  and  negative  mass  Instabilities  for  this  design,  concluding  that 
the  beam  breakup  Instability  poses  few  problems  for  Induction  module  0's  of 
order  six,  a  realistic  value.  The  negative  mass  Instability  was  found  to 
be  somewhat  more  threatening,  although  there  was  reason  to  hope  that 
thermal  effects  would  reduce  the  predicted  growth  rates,  especially  for 
large  toroidal  mode  numebrs.7  Our  work  was  based  on  the  dispersion  rela¬ 
tion  of  Ref.  2.  Here,  we  repeat  the  analysis  using  Eqs.  (23)  -  (27). 

We  treat  the  racetrack  cavity  as  a  torus  of  major  radius  R  =  70  cm. 
(Straight  sections  of  the  racetrack  are  expected  to  have  a  favorable,  but 
very  small,  effect  on  total  negative  mass  instability  growth.)  Varying  the 
assumed  radius  has  little  affect,  because  changes  In  the  growth  rate  are 
approximately  balanced  by  changes  In  the  path  length.  In  accordance  with 
the  accelerator  design,  the  cavity  and  beam  minor  radii  are  taken  to  be  b  = 
7  cm  and  a  *  1  cm.  A  Be°  *  2  kG  (u^  *  1.73)  guide  field  is  applied. 

The  numerically  determined  growth  rate  for  these  parameters  and  i  = 

13,  a  toroidal  mode  number  near  the  upper  end  of  the  range  for  which  the 
model  is  valid.  Is  plotted  as  the  solid  curve  In  Fig.  2.  The  result  from 
Ref.  6  is  shown  as  a  dashed  line.  The  average  growth  rate  Is  reduced  by 
about  a  factor  of  two.  Our  more  optimistic  findings  suggest  a  total  nega¬ 
tive  mass  Instability  amplification  factor  of  fifteen  e-foldlngs.  This 
growth  level  probably  can  be  reduced  to  an  acceptable  level  (say,  five  e- 
fol dings)  by  thermal  effects  or  design  changes. 
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NEW  RESULT 


Figure  2.  Predicted  negative  mass  instabi 
of  energy  for  the  ONR  racetrack 
with  a  2  kG  guide  f ield. 


Figure  3  Illustrates  the  corresponding  Be°  *  0  growth  rates.  Peak 
values  are  reduced  somewhat  relative  to  the  earlier  results,  but  the  high 
energy  cutoff  Is  shifted  up  in  energy  by  a  factor  of  twenty.  For  y  >  20, 
Eq.  (34)  reproduces  the  solid  curve  well. 
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1.  INTRODUCTION 


The  modified  betatron  concept,1-3  illustrated  in  Fig.  1,  may 
provide  a  compact  means  of  accelerating  Intense  electron  beams  to  high 
energies.  A  dispersion  relation  for  the  linear  stability  of  the  electron 
ring  In  the  device  has  been  derived  by  Sprangle  and  Vomvoridis.14  In  this 
report,  we  show  that  some  of  the  approximations  In  their  derivation  are  not 
well  justified,  and  we  obtain  more  accurate  expressions.  In  Sec.  II,  the 
approximation  that  the  phase  velocity  of  unstable  waves  is  approximately 
the  same  as  the  beam  velocity,4  V^*Vb,  Is  discarded.  This  signifi¬ 
cantly  alters  the  results  obtained  In  two  ways.  Firstly,  the  growth  rates 
obtained  are  typically  two  to  ten  times  larger.  Secondly,  we  find  that  the 
conventional  negative  mass  Instability  does  not  exist  in  modified  beta¬ 
trons.  Rather,  the  beam  Is  subject  to  a  predominantly  transverse  instabil¬ 
ity  at  high  energies.  Me  have  made  a  rough  estimate  of  the  effect  of  a 
spread  In  beam  energy  on  this  mode.  In  Sec.  Ill,  we  examine  the  effect  of 
a  moderate  spread  In  beam  energy  on  the  transverse  resistive  wall  Instabil¬ 
ity.  We  find  that  In  some  cases,  the  effect  Is  negligible  because  Vb  - 
Is  too  large.  For  the  most  dangerous  nonresi stive  and  resistive 
instabilities,  however,  significant  damping  Is  expected. 
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II.  LINEAR  DISPERSION  RELATION 


A.  Derivation 

Our  analysis  follows  that  of  Ref.  4,  except  that  we  assume  a 
monoenergetl c  beam.  The  details  of  the  derivation  are  given  in  the  Append¬ 
ix,  and  here  we  give  only  the  main  points.  The  beam  is  modeled  as  a  circu¬ 
lating  ring  of  charge  which  can  displace  rigidly  in  the  transverse  direct¬ 
ion  and  which  can  compress  In  the  toroidal  direction  (see  Fig.  1). 


In  equilibrium,  the  beam  is  positioned  at  the  center  of  the 
minor  cross-section  of  the  torus,  and  executes  a  cyclotron  orbit  in  the 
mirror  Bz  field.  Toroidal  corrections  to  the  field  equations  are 
dropped,  so  that  the  m  *  0  and  m  *  1  fields  are  not  directly  coupled.  They 
are,  however,  coupled  via  the  perturbed  charge  and  current.  Thus,  the  m  = 

0  component  of  the  charge  density  p  satisfies 


where  r(e)  Is  the  radial  location  of  the  center  of  the  beam,  and  Vr,  Ve 
are  the  beam  velocity  components.  The  second  term  in  Eq.  (1)  shows  that  a 
rigid  transverse  (m  *  1)  displacement  contributes  to  the  perturbed  net 
(m  *  0)  charge  density.  Contributions  from  perturbed  m  =  0  quantities  to 
the  m  *  1  charge  density  are  second  order  In  the  beam  transverse  displace¬ 
ment,  and  so  do  enter  the  linear  dispersion  relation.  Consequently,  the 
perturbed  m  «  1  fields  can  be  computed  directly  in  terms  of  the  transverse 
displacements  of  the  beam.  The  results  are  substituted  Into  the  m  *  0 
field  equation  for  the  perturbed  toroidal  electric  field  Eg!1),  namely 
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(see  Appendix  for  definitions  and  normalizations.)  Linearizing  Eq.  (1),  we 
oh*’’r- 


Solving  this  equation  with  appropriate  conducting  wall  boundary  conditions 
yields  the  linear  dispersion  relation 
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(1  +  2in  a/rb)  (1  -  (1  +  1)  E|i) 


This  equation  differs  from  the  results  of  all  earlier  work  in  that  the 
approximation  u>  -  iflZ0/y0  has  not  been  roade*  Also,  the  first  term 
on  the  second  line  is  new. 


B.  Nonresistive  Instabilities 


Equation  (3)  has  some  unstable  roots  due  to  the  coupling  of 
longitudinal  and  transverse  modes  of  oscillation.  The  instabilities  per¬ 
sist  when  the  wall  conductivity  is  infinite.  The  instabilities  are  low 
frequency  in  the  sense  the  transverse  component  of  their  motion  is  associ¬ 
ated  with  the  slow  rotation  frequency  wg  «  “rta>z/(n0O/^o) • 

The  u.am  can  also  oscillate  transversely  at  the  fast  rotation  frequency, 
06o/yo»  bub  there  are  no  nonresistive  instabilities  associated  with 
this  resonance. 
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We  can  clarify  the  origin  of  the  nonresl stive  Instabilities  by 
simplifying  Eq.  (3).  We  assume  <  flg0/Yo*  Aw2  <  w*  and 
obtain 


p  5  J  porb  (1  +  2ln  a/rbJ 
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where  a  *  yo/U  -  n  -  ns  rb*/aO.  The  function  P(Aoj)  has  a  dif¬ 
ferent  character  depending  on  whether  |a>Bj  <  *ftzo/Yo  or  |wbI 
>  to  /yo»  »s  shown  In  Fig.  2.  For  typical  betatron  parameters,  the 
point  juB|  *  tfl20/Yo  occurs  approximately  at  y0  =  Ytran  * 

[4v  r02/ a2]1/3  where  v  is  Budker's  parameter.  When  yQ  <  Ytran, 
the  roots  of  the  quartic  P(auj)  =  1  consist  of  two  complex  conjugate  pairs. 
For  y0  >  y tran  ’  we  have  tw0  rea^  roots  and  a  complex  conjugate 
pair.  We  note  that  the  conventional  negative  mass  Instability5  is  not  pre¬ 
sent  In  typical  modified  betatrons.  The  derivation  of  the  dispersion  rela¬ 
tion  for  the  latter  Instability  involves  the  replacement  of  y^2  *  (1  - 
«2r02/i2)-1  by  yb2  *  (1  -  V^2)-1  In  the  field  equation.  This  proced¬ 
ure  Is  valid  only  if  2vy0  (1  +  2m  a/r 5)  <  1  which  Is  not  the  case  for 
modified  betatron  parameters  (cf.  Table  1). 


Frequencies  and  growth  rates  for  the  1  -  1,2,  3,  and  4  nonre- 
sistive  modes  obtained  by  solving  Eq.  (3)  numerically  for  the  parameters  in 
Table  1  are  given  In  Figs.  3,  4,  and  5.  As  we  have  seen,  for  y0  < 
ytran  there  are  tw0  unstable  modes,  one  with  u  >  *ft20/ Yo»  the 
other  with  u  <  tnzQ/y0,  and  we  term  these  modes  "fast"  and  "slow" 
accordingly.  In  Figs.  3  and  4,  only  the  slow  modes  are  depicted  for  clar¬ 
ity.  The  maximum  growth  rates  are  for  y0  >  Ytran.  and  since  most  of 
the  acceleration  period  lies  in  this  region,  we  shall  examine  the  region 
more  closely.  For  a  «  2y02  >  1,  Eq.  (3a)  reduces  to 

Au  (Au>2  -  «g2)  +  v  (1  +  2in  a/r^)  wau>g2/Y0^  *  0  (3b) 
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The  condition  for  this  cubic  in  tu  to  have  complex  roots  is  v(l  + 

2tna/rjj)  uo/y03  >  2  «b/(3/3).  This  criterion  yields  the  upper  bound 
on  the  unstable  range  of  yo*  namely 

Ymax  '  (1  ♦  2m  ./r„)  « 

where  we  have  assumed  wg  <  in20/Yo*  This  expression  gives  rmax 
*  153  for  the  parameters  in  Table  1(a).  The  exact  numerical  results  give 

w  -  156- 

For  y02  <  Ymax2  the  complex  roots  of  Eq.  (3b)  are  given 
approximately  by 

A«  -  Cy0v  1(1  +  2  in  a/rb)/(2r^B^)31/3  e1*  , 

where  $  =  2n/3,  4ir/3.  Thus,  the  growth  rate  scales  as  t1/3,  B~2/3, 
etc.  For  the  parameters  in  Table  1(a),  this  expression  yields  w  =6.81xl0"3 
+  i  2.4x10'**  for  70  *  50,  compared  to  the  exact  answer  u>  =  6.82xl0*3  + 

1  2.1x10'**.  Numerically  we  find  that  throughout  most  of  the  range  of  this 
instability  we  have  *  wg,  so  that  the  mode  is  mostly  transverse  in 
character.  The  conventional  negative  mass  instability  is  longitudinal  in 
character,  being  associated  with  the  Lw  -  0  resonance. 

A  comparison  between  our  dispersion  relation,  Eq.  (3)  and  that 
in  Ref.  4  is  given  in  Fig.  6.  The  mathematical  differences  between  the  two 
dispersion  relations  were  described  in  Sec.  1 1  A.  Equation  (3)  gives  growth 
rates  which  are  two  to  ten  times  larger  than  those  from  Ref.  4.  We  discuss 
che  effect  of  a  thermal  spread  in  energy  on  these  instabilities  in  subsec¬ 
tion  D  below. 

C.  Resistive  Wall  Instabilities 

The  presence  of  resistive  material  in  the  walls  of  the  betatron 
gives  rise  to  additional  Instabilities,**  and  modifies  the  growth  rates  of 
nonreslstive  Instabilities.  To  Illustrate  this  effect,  we  have  chosen  a 
stainless  steel  wall,  for  which  the  conductivity  a  is  5.2xl06  in  normalized 
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units  (see  Appendix).  The  results  for  the  %  *  1  mode  are  shown  in  Fig.  7. 
The  resistive  wall  has  little  effect  on  the  nonresistive  instabilities. 
However,  some  modes  whose  growth  rates  are  zero  for  0  =  •  are  driven 
unstable  by  the  resistivity.  They  are,  the  fast  mode  in  the  region 

Y0  <  >tran  an<*  a  s1ow  mode  re9^on  To  >  I'tran’  denoted 

by  A  and  B  respectively  in  Fig.  7(a).  (In  Fig.  7(b),  branch  A1  is  unstable 

even  for  g  -  ■>.)  Branch  B  Is  due  mainly  to  the  term  en  in  Eq.  (3).  The 
growth  rates  of  this  branch  are  much  smaller  than  those  obtained  by  using 
the  approximate  en  in  Ref.  4.  Since  the  resistive  modes  are  driven  by  the 
boundary  condition  at  the  wall,  they  are  sensitive  to  the  value  of  a,  the 
minor  radius  of  the  torus.6  This  is  why  the  resistive  mode  growth  rate  is 
smaller  in  Fig.  7  than  in  Fig.  6  (cf.  Table  1).  The  growth  rate  is  approx¬ 
imately  Independent  of  y  . 

The  slow  mode  associated  with  the  toroidal  magnetic  field  cyclo¬ 
tron  resonance,  u  -  Is  also  driven 

ZU  0  00  0 

unstable  by  wall  resistivity.6  As  Indicated  in  subsection  B,  none  of  the 
modes  associated  with  this  resonance  are  unstable  when  0  s  ®.  With  finite 
wall  conductivity  the  mode,  which  is  primarily  a  transverse  oscillation, 
becomes  unstable  when  u  goes  through  zero  and  becomes  positive.  In  a  beta¬ 
tron,  ozQ  «  y0  during  the  acceleration,  so  that  the  instability 

turns  on  when  n  „  *  £!„„/£  and  continues  for  the  remainder  of  the 
zo  eo 

acceleration  period.  This  behavior  Is  shown  in  Fig.  8.  Again,  the  differ¬ 
ence  in  growth  rates  between  the  two  parts  of  the  figure  is  due  mainly  to 
the  differences  In  the  quantities  a  and  B._  in  Table  1.  For  large 
Y0,  the  growth  rate  Is  approximately  Independent  of  yq. 

D.  Practical  Implications  for  Betatrons 

For  the  sample  parameters  given  in  Table  1,  it  is  clear  that  the 
nonresistive  Instability  In  the  region  yo  >  Ttran  is  the  most  import¬ 
ant  Instability.  Thus,  for  the  parameters  In  Table  1(a),  the  number  of  e- 
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foldings  of  the  i=l  component  during  a  1  millisecond  acceleration  time  is 
about  4000.  This  result  is  for  a  monoenergetic  beam,  and  gives  an  upper 
bound  on  the  growth.  We  can  estimate  the  effect  of  a  spread  in  beam  energy 
as  follows.  The  thermal  spread  enters  the  model  in  the  combination  w  - 
1(oZo/t°  -  kAP0),  where  aP0  is  the  spread  in  canonical  toroidal 
momentum  (cf.  Eq.  (4)).  The  instability  for  y0  >  Ytran  is  associated 
with  resonance  tu  =  wg.  Therefore,  a  small -thermal -expansion  for 
this  mode  is  an  expansion  in  the  parameter  e2=  UkAP0)2/(Au>  -  wg)2. 

If  c2  <  1,  Landau  damping  is  negligible,  whereas  if  e2  ^  1,  we  expect  sig¬ 
nificant  damping.  As  an  example,  we  use  the  numerical  results  shown  in 
Fig.  3(a),  and  assume  an  initial  spread  in  y0  of  5%.  Then,  for  the  t=l 
mode  at  y0  =  50,  we  obtain  e2  »  4,  so  that  we  can  expect  a  significant 
reduction  in  growth  rate.  A  more  rigorous  treatment  of  thermal  effects  is 
needed  to  confirm  this  result. 
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III.  THERMAL  EFFECTS  ON  RESISTIVE  INSTABILITIES 

It  has  been  suggested  that  a  moderate  spread  In  beam  particle 
energies  may  reduce  Instability  growth  rates  to  acceptably  low  values 
through  Landau  damping.  Here  we  look  at  the  effect  of  a  thermal  spread  on 
the  transverse  cyclotron  resistive  wall  instability.  We  choose  this  case 
because  the  dispersion  relation,  Eq.  (4)  is  relatively  simple  and  does  not 
require  numerical  solution. 

A.  High  Frequency  Limit:  6/a  <<  1 

From  Ref.  4  the  approximate  dispersion  relation  for  the  cyclo¬ 
tron  mode  including  thermal  effects  is 

1  *  “s/  <«' 

“r  "  Aa)  1  ^BO/yO 
where,  in  normalized  units. 


Au>  a  a)  -  t(n20/y0  -  kaP)  , 


2  1  2  2 

wr  (•£  -  n$)  nz0/yQ  .  (n  *  1/2  is  assumed) 


g(aP)  *  distribution  function  of  toroidal  canonical  momentum 
spread.  See  Appendix  for  additional  definitions. 
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Write  Aw2  -  wr2  ?  6ti>nQ0/y0  *  (Aw  +  IkAP  -  ox )  (aw  +  £kaP  -  a2).  Assuming 
|o2|»lai],  the  Instability  comes  from  the  following  choice  of  roots. 


w_  Y- 
r  o 


=  ‘  neo/Yo 


For  g,  choose  a  flat-topped  distribution  function, 
gUP)  -  for  |aP|  <  aPq  , 


g(AP)  *  0  for  j aP j  >  aP£ 


Performing  the  integration  in  Eq.  (4),  we  get 


2£kAPQ(a1  -  a2) 


AwQ  -  tkAP0  -  a2 
Acj_  +  £kAP  -  bo 


Aa)0  +  £kAP0  -  oj 
Au,0  ‘  UAP0  -  oj 


where  Aw.  s  w  -  in  J y*  The  mode  we  are  concerned  with  has 
o  zo  o 

Aw  «  02 .  To  do  a  small -thermal -spread  expansion,  we  assume  kAPQ<<  Aw  -  a2 

In  what  follows,  we  shall  in  essence  be  checking  the  consistency  of  these 
two  approximations.  Expanding  Eq.  (7),  we  obtain 


i --4  ( r *  1 


(£kAP0)C 
(Awq  -  a2)' 
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Assuming  6/a  «  1,  the  real  part  of  the  frequency,  awr.  Is 


Aur  *  '  fl8o/Y 


0 


y  (Ao>r-  ag )  J 


Our  expansion  parameter  Is  thus 


Au»r  -  <*2 


o  eo  'o 


With  aPq  =  Ythr0.  and  yQ  large  enough  such  that  n$  <<  1,  we  have 

4i*Wo  Yth 


If  e  <<  1  for  a  given  choice  of  parameters,  then  our  small- 
thermal  -spread  expansion  Is  valid.  In  this  limit,  there  is  no  Landau  damp¬ 
ing  from  a  flat-topped  distribution.  If  e  >>  1,  on  the  other  hand,  the 
phase  velocity  of  the  mode  lies  well  within  the  distribution  of  particle 
velocities.  The  mode  Is  then  highly  damped.  Putting  in  numbers  from  Table 
l(a)with  yQ  =  50,  we  obtain  e  =  3.3*103 (t^^/t0) •  Assuming  a 
10%  spread  in  at  the  beginning  of  the  acceleration  period,  we  have 
Yth  *  Thus  c  *  33  at  Y0  *  5®*  Consequently,  there  will  be 

significant  Landau  damping  of  this  mode. 

B.  Low  Frequency  Limit:  u  »  0 

When  u  *  0,  the  small-thermal-spread  expansion  of  Eq.  (7)  leads 


3/2  “so  in/4  1  ,  n 

u  +  -  n  ®  V 1  +  t  e  J  =  u 

a£  A 


wheren^  *  n$(l  -  r2/a2)  ,  and  n  =  B^(r2/a3)(2/o)1/Z.  The 


10 


unstable  solutions  to  Eq.  (9)  are 


„  2 
“to  " 

°ec/To 


h 


where  ,  *  n/6,  -  7,/6.  In  this  case,  we  obtain 


7th  ta2  /“ ao\2/3  /oV 

?T  ^  W  W 


For  the  parameters  in  Table  1(a),  we  obtain  e  »  27 ( .  Thus, 
for  an  initial  10%  spread  in  y  there  will  be  substantially  less 
Landau  damping  in  this  case  than  where  6/a  <<  1.  For  the  parameters  in 
Table  1(b),  e  *  case*  Landau  damping  will  be 

negligible,  and  Eq.  (10)  shows  that  there  will  be  a  slight  increase  in  the 
growth  rate  due  to  thermal  effects. 


IV .  SUMMARY 


We  have  rederived  the  dispersion  relation  for  linear  instabil¬ 
ities  in  betatrons  based  on  the  simple  model  of  Ref.  4.  Our  analysis  shows 
that  at  high  beam  energies,  the  dispersion  relation  does  not  reduce  to  that 
of  the  conventional  negative  mass  Instability.  Instead,  we  find  a  mode 
which  is  primarily  transverse  in  nature.  Furthermore,  we  obtain  growth 
rates  which  are  from  two  to  ten  times  larger  than  those  obtained  in  pre¬ 
vious  calculations.  We  have  estimated  the  effects  of  a  moderate  spread  in 
beam  energy  on  nonresistive  and  resistive  instabilities,  and  find  that  sig¬ 
nificant  damping  Is  expected.  A  more  rigorous  calculation  is  needed  to 
prove  this. 
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APPENDIX 

DERIVATION  OF  DISPERSION  RELATION 

Figure  1  illustrates  the  physical  parameters  of  the  system.  The 
beam  is  modeled  as  a  line  of  charge.  In  equilibrium,  it  is  situated  at 
r  *  r0  2  *  0,  and  executes  a  cyclotron  orbit  in  the  mirror  Bz  field. 

As  pointed  out  in  Sec.  IIA,  only  the  transverse  motion  of  the  beam  needs 
to  be  considered  when  computing  the  perturbed  transverse  fields.  When  the 
center  of  the  beam  is  displaced  rigidly  to  position  {r0  +  tr,  lz) ,  it 
experiences  the  following  fields:4 


Applied  fields:  Bz  =  B2Q  (1  -  nar/rQ)  , 

Br  *  -®20  "lr/r0  ■ 


8e  ■  Beo  11  -  4r/ro>  ■ 


Induced  fields:  E. 


"o  rb 

A  r  , 
a 


•  d)  .  no  rb 


T  “2  *z 

a 


B'11  -  ^  rb.. 


-  (1  -  r2/a2)  *  Ur,  (A2) 


where 


• a  4 11  ttY  (i?) 
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In  our  normalization  scheme,  frequencies  are  normalized  to  which  is 
defined  by  c/u>  *  1  cm.  Lengths  are  normalized  to  c/w  ,  velocities 

®  o  p  ° 

to  c,  fields  to  mo*)0/e,  and  densities  to  w£m/4ire  ,  where  m  and  e 
are  the  electronic  mass  and  charge  respectively.  Thus,  for  example, 

B20  and  flz0  have  the  same  normalized  values.  The  conductivity  is 
normalized  tou>0/4ir.  In  Eqs.  (Al)  and  (A2)  above  n  is  the  external 
field  Index,  l.e.,  rQ  3/3r  in  Bz(r0),  n0  Is  the  equilibrium  beam 
density  and  V  Is  the  equilibrium  azimuthal  beam  velocity.  A  positive 


eo 

beam  charge  Is  assumed. 

A  . 

•  y  -  + 

T  r 


The  equations  of  motion  of  a  beam  particle  are. 


dpr 

Si 


E  +  V  B  -  V  B 
r  ez  2  6 


dpz 

■at 


E  +  V  B  -  V  B 
z  re  e  r 


V  V 
e  r  . 
y  ■—  + 


dp 


_e 

dt 


E  +  V  B  -  V  B 
e  z  r  r  z 


where  (Vr,  Ve ,  Vz)  and  (pr,  pg,  pz)  are  the  velocity  and  momen¬ 
tum  components  of  a  beam  particle.  Linearizing  these  equations,  we  obtain 

lt  *  <4  -  Thr  .  ii  -  -  4  ^  , 

'  To  ’0 


az  +  (w*  -  0  az  -  af 

*  '  o 


0  , 


(A3) 


<JV 


(I) 
e _ 

sr~ 


(1) 
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where  the  superscript^  denotes  perturbed  quantities, 

\  rb/a2>eL/r0'  “I  *  <"  -  ",  rb/iJ)Bfo/',0-  "s  ■  "o/(2l(0B2ol 
and  e  •  nsT2v|0(l  -  r2/a2)5B2Q/'f2.  Assuming  ar,  u  -  exp  {-iut  ♦  lee), 
we  obtain  the  following  solutions  to  Eqs.  (A3), 


2  (aw2  -  w2  ♦  5)  B  h 

e  To  '  z  zo' To 

,  - B - 

ve  )yo  ^Bzo^Yo^®eo^Yo^ 

"  5 


p  2—  2  2—  222 

where  Aw  *  u  -  *vet)/r0  and  D  =  (Aw  -  wj;  +  c)(Aw  -  w^+  0  -  Aw 
Using  the  third  member  of  Eqs.  (A3),  V.U)  can  be  expressed  in  terms 

O 

of  E_ ( 1 ) ,  the  only  unknown.  To  close  the  system,  we  obtain  a  field 

9 

equation  for  Ee ( 1 ) .  We  assume  that  the  beam  excites  only  the  m  =  0 
component  of  E  (*).  This  is  reasonable  provided  rj,  «  a,  since 
higher  m  number  components  go  to  zero  at  the  center  of  the  minor  cross 
section  of  the  torus.  Further,  we  assume  that  only  the  lowest  radial  mode 
Is  excited,  so  that  the  eigenfunction  Is  approximately  constant  over  the 
beam  cross  section.  This  is  valid  provided  jt/r0j  -  jwj  <<  2n/a.  Then 
from  Maxwell's  equations  we  obtain 


V2  E*1* 

\L  Le 


Up 


(A5) 


where  "1"  refers  to  the  transverse  direction,  is  the  perturbed 
charge  density,  and  i}\  s  the  perturbed  azimuthal  current.  Since  the 
eigenfunction  Is  assumed  to  be  flat  In  the  center,  the  perturbed  charge 
density  p^  Is  proportional  to  the  perturbed  line  charge  v^.  To 
obtain  an  expression  for  we  use  the  continuity  equation  for  p. 
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FT? u  v A J.'UI 


: 


v 


* » 


ip/it  +  V  •  (pV)  *  0. 


(A6) 


Put  p  »  v€(R  -  R0(e))  6{*  -  ♦0(®^/R»  where  R»  ♦  are  local  cylindrical 

coordinates  (r  -  r0  «=  Rcos$,  z  *  rs1n*>),  and  (R0,  *Q)  is  the 

position  of  the  displaced  beam.  Multiplying  Eq.  (A6)  by  JRdRd$,  we  obtain 


vV  vV 

lv  +  !_r  +  L.  (Li)  =  o 

at  r  de  r 


(A7) 


LV 

[v“ 

LV 

LV 

LV 


Linearizing,  and  replacing  v  by  p,  we  have 


111  *  <t»ove1>/ro  ‘  »o"rlU/ro’/4“ 


^  ■  .«VM‘  *  r%0„< 


Sc 

u 


where  p  is  the  unperturbed  charge  density  .  Equation  (A5)  becomes 


o2f(1)  - 
vJ.Le 


,  F{1) 

*poEe 
— TT 

YAo.  rQ 


i 

~r 

yn 


uB  /v 
zo'  To 


(a)2  -  Au>2  -  O 


1  - 


2  2 
“  ro 
~~T 


IN 


“  poBzo/yo  #  2 
to  yam  Uz 


Ato  - 


J 


£  AE 


(1) 


(A8) 


* 
*  « 


To  obtain  the  dispersion  relation,  we  need  to  solve 

(1) 


72  E(1) 
yl  Le 


AE 


r  <  r 


72E(1) 

viS 


r.  <  r  <  a 
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”iE“’ 


-iowE 


(1) 


r  >  a, 
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TABLE  1.  SAMPLE  PARAMETERS  FOR  MODIFIED  BETATRON 

The  values  for  B.  are  approximate  practical  upper  and  lower  bounds. 
00 


IUANTITY  SYMBOL  VALUES 


(a) 

(b) 

Major  Radius 

ro 

150  cm 

100  cm 

Minor  Radius 

a 

5  cm 

10  cm 

Beam  Radius 

% 

1  cm 

1  cm 

Toroidal  Magnetic  Field 

Beo 

6  (10  kG) 

1.5  (2.5  kG) 

Beam  Current 

V 

0.59 

0.59  (10  kA) 

Beam  y-factor 

7o 

5-100 

Transition  Energy 

ytran 

12.9 

6.2 

Acceleration  Time 

Ta 

3  x  107  cm 

(1  millisec) 
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Figure  1.  Illustration  of  modified  betatron  concept.  The  major  radius  of 
the  torus  Is  r0,  the  minor  radius  is  a  and  the  beam  radius  is 
tk.  The  external  magnetic  fields  consist  of  a  focusing  mirror 
field  B2  and  a  toroidal  field  Be. 


20  0 


Figure  2.  Illustration  of  nonreslstlve  Instabilities  in  betatrons. 

The  dispersion  relation  Is  P»l.  The  roots  are  denoted  by 
Tj  ,  r2,  etc.,  and  brackets  (,)  denote  complex  conjugate 
pairs.  In  (a)  and  (b),  we  depict  two  regimes  of  Instability 
In  the  modified  betatron.  In  (c),  we  show  the  origin  of  the 


longitudinal  negative  mass  Instability,  which  requires 
k  •  2vyQ  (1  ♦  2  tn  a/rb)  <  1. 


r 
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NL 


Figure  3.  Real  and  Imaginary  parts  of  the  nonresi stive  l  •  1,  2,  3,  4 
inodes  obtained  by  solving  Eq.  (3)  numerically  for  the  equil- 

in  Table  1(a).  The  frequencies  are 
»  <  Ytran •  only  the  slow 


Ibrlum  parameters  given  ir 
In  units  of  3x10*°  sec*lt  For 
is  Included  (cf.  Fig.  2). 
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Growth 


betatron  parameters 


Figure  6.  Comparison  between  growth  rates  obtained  from  Eq.  (3)  (Curve  A) 
and  those  obtained  from  the  dispersion  relation  In  Ref.  4  (Curve 
B),  for  the  nonreslstlve  t  -  1  Instability.  The  betatron  param¬ 
eters  are  from  Table  1(a). 
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Figure  7. 


Growth  retes  for  the  i  »  1  Instability  with  perfectly  conducting 
walls  (dashed  lines)  and  stainless  steel  walls  (solid  lines). 

Part  (a)  Is  for  the  parameters  In  Table  1(a),  and  part  (b)  Is  for 
parameters  In  Table  1(b).  Branches  A  and  B  are  modes  which  have 
become  unstable  due  to  the  wall  resistivity  alone.  In  part  (b) 
the  solid  and  dashed  lines  are  Indistinguishable  (the  growth  rate 
of  Branch  B  Is  approximately  4x10-6  cm-  ) .  Branch  A*  Is  unstable 
for  even  o  ■ 


ft  growth  rates  of  the  transverse  resistive  wall  cyclotron  mode. 
"• "*  *•  *?!(,)  "for  th.  in  Tnblt  !(•)  nnd  port  b)  Is 

for  thoit  in  Tibi*  1(b).  for  (.),  I  ■  »  onb  for  (b)  t  •  S. 

The  instability  turns  on  when  R*o ■  k  is 

(•  45  for  case  (a)).  The  height  of  the  Initial  peak  is 

Independent  of  t  and  to  (c* •  **f.  6'* 
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Brendan  B.  Godfrey  and  Thomas  P.  Hughes 
MISSION  RESEARCH  CORPORATION 


ABSTRACT 

Resistive  wall  Instabilities  In  modified  betatrons  are  analyzed  In 
several  limits.  The  moderate  frequency,  negative  energy,  m*l  spacecharge 
and  cyclotron  waves  are  found  to  be  most  dangerous,  potentially  capable  of 
disrupting  acceleration  for  typical  betatron  parameters.  A  moderate  spread 
In  electron  energy  can,  however,  stabilize  these  modes. 


Spacecharge  limits  the  circulating  current  In  a  conventional  betatron 
to  tens  of  amps,  unless  a  very  high  Initial  electron  energy  Is  assumed. 
Adding  a  strong  toroidal  magnetic  field,  however.  Increases  this  limit  by 
some  two  orders  of  magnitude. The  toroidal  field  also  Improves  beam 
stability,4  while  a  higher  current  has  the  opposite  effect.  This 
report  discusses  resistive  wall  Instabilities5  In  a  high  current  elec¬ 
tron  betatron  with  an  applied  toroidal  magnetic  field,  referred  to  as  a 
modified  betatron. 

Resistive-wall -driven  Instability  occurs  for  both  slow  spacecharge  and 
slow  cyclotron  beam  modes,  and  over  a  wide  range  of  frequencies.6 
Approximate  analytical  growth  rate  formulas  are  obtained  In  this  report  for 
the  various  branches  of  the  resistive  wall  Instability  on  a  cold  beam.  The 
growth  rate  expressions  are  then  evaluated  for  a  typical  set  of  modified 
betatron  parameters.  He  find  the  spacecharge  and  cyclotron  wave  Instabil¬ 
ities  at  frequencies  comparable  to  the  electron  circulation  frequency  to  be 
the  most  dangerous.  Although  slowly  growing,  they  persist  through  a  large 
fraction  of  the  acceleration.  If  one  can  produce  a  spread  of  a  few  percent 
In  electron  energy,  however,  the  modes  can  be  stabilized.  Alternatively, 
they  can  be  avoided  by  limiting  acceleration  times  to  tens  of  micro¬ 
seconds. 

The  structure  of  the  resistive  wall  Instability  dispersion  relation  Is 
obtained  for  highly  conducting  cavity  walls  by  a  perturbation  expansion 
(formally,  a  Raylelgh-Rltz  approximation)  about  the  beam  modes  In  a  per¬ 
fectly  conducting  cavity.  This  procedure  yields 

2(«-«0)U  *  P, 


(1) 


0 


where  u0  Is  the  unperturbed  wave  frequency,  U  Is  the  wave  energy,  and  P 
Is  the  outward  Poyntlng  flux  at  the  cavity  wall  due  to  finite  conductiv¬ 
ity.  P  Is  evaluated  easily  In  the  long  wavelength  limit  to  be7 

P  *  RuSB2  eU/4  .  (2) 

Here,  R  Is  the  minor  radius  of  the  toroidal  cavity,  B  is  the  wave  magnetic 
field  at  the  cavity  wall,  and  5  is  the  skin  depth, 

6  *  /z/au  •  (3 ) 

Frequencies  jnd  conductivities  are  expressed  in  inverse  cm;  the  speed  of 
light  is  unity.  Also,  a  factor  of  4*  is  absorbed  into  the  conductivity  c. 

Solving  Eq.  (l)-(3)  for  the  perturbed  frequency  yields 

u)1/2(a.  -  u>0)  -  2(U/u,)"1(u1/26)RB2  eW4  .  (4) 

The  right  side  of  Eq.  (4)  Is  approximately  independent  of  frequency, 

O 

because  U  «  u  for  beam  modes.  (This  assertion  is  incorrect  for  m  =  0 
spacecharge  waves.  Nonetheless,  the  resulting  dispersion  relation  has  been 
found  nunerically  to  apply  reasonably  well  even  then.7)  Evaluating  the 
wave  energy  to  complete  the  derivation  is  straightforward  but  tedious. 
Instead,  we  extract  the  needed  term  from  the  recent  work  of  Sprangle  and 
Vomvoridis.9 


be 


These  authors  determined  the  Imaginary  part  of  u  (the  growth  rate)  to 
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n$yfia2/R3 
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(5) 
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with 


A2  -  1  ♦  2(1  -  2nsa2/R2)Qz2/ne2  (6) 

for  Moderate  frequency  m»l  waves.  The  beam  radius  is  a.  Otherwise,  sym¬ 
bols  are  defined  as  In  Ref.  9.  Recasting  Eq.  (5)  as 

2.  2  /p3 
ea  /R 


we  see  that  toroidal  corrections  to  the  cylindrical  drift  tube  dispersion 
relation  enter  only  through  a.  Comparison  of  Eqs.  (4)  and  (7)  yields 


1/2/  x 

U  lb)  “  Up) 


u2.2/r3 


a  /  . 

0  fieA 


For  Increased  generality  we  have  included  a  geometrical  factor  f  by 
analogy  with  cylindrical  drift  tube  results.7 


f  -  (HI2  +  t2R2/R  2 )  r, - 

1  0  I l*  +  1 


«. ,  1  “"V 
W„1  +  V-'  l"1 


R0  Is  the  major  radius  of  the  torus,  I  and  K  are  modified  Bessel 
functions,  and  t  amd  m  are  the  toroidal  and  poloidal  wave  numbers.  With 
a/R  and  tR/Rn  both  small,  f  Is  aproxlmately  1  for  m  *  1  and  varies  as 
m2(a/R)2  m_1*  for  m  >  1.  The  Instability  Is  weaker  for  m  =  0  where 
f  falls  off  as  (*a/R0),2  and  for  m  <  0  (wave  poloidal  hell  city  opposite 
that  of  the  electrons),  where  f  decreases  as  m(a/R)  ' 


For  slow  cyclotron  waves,  characterized  by 
u  •  i  nz/y  -  n0/Y  . 


the  peak  growth  rate  from  Eq.  (8)  occurs  at  w0  »  0. 


2„2/»3 

_  -1/3  “p  a  /R  f 

“  °  [  v  *J 


The  constraint  u>q  «  0  Implies  t  »  Qe/fl2.  During  acceleration,  the 
vertic-  magnetic  field  nz  (and  the  beam  energy  y,  which  is  proportional 
to  n2)  increases  adiabatically.  Hence  u0  •  0  is  satisfied  for  only  a 
limited  time  period. 


?!  i*L 
v  Ut 


Equating  the  instability  band  width  Au>  to  the  growth  rate  r,  replacing 
by  in2,  and  solving  for  at,  we  obtain  the  total  growth  occurring  as  the 
instability  passes  through  resonance 


rz  v 

r  At  "  tlaz/ y)  dy/dt 


Eq.  (13)  :s  maximized  late  In  the  acceleration  cycle,  when  y  is  large  and 
the  resonant  t  value  Is  small. 

After  passing  through  u0  -  0,  the  cyclotron  resistive  wall  instabil¬ 
ity  for  flxedt  does  not,  of  course,  cease  but  Instead  transitions  smoothly 
to  a  lower  growth  rate  regime. 


.  -1/2  -1/2 
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Setting  m  *  1  recovers  Eq.  (5).  Although  the  growth  rate  of  Eq.  (14)  is 
somewhat  smaller  than  that  of  Eq.  (11),  the  corresponding  total  growth  of 
the  former. 


(rat)eff  -  (2o) 
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maty  exceed  Eq.  (13).  The  Integration  In  Eq.  (15)  extends  from  the  time  at 
which  «0  Is  greater  than  a  few  times  the  u  In  Eq.  (11)  until  the  end  of 
the  acceleration. 


Parallel  Instabilities  occur  for  slow  spacecharge  waves,  characterized 


wo  *  ~  '  (16 ) 

where  Is  the  beam  pololdal  rotation  frequency,  a  small  quan¬ 
tity.®  For  growth  to  occur  at  the  rate  In  Eq.  (11),  we  must  have  m  - 
tflz/ywB*  Even  1  as  sma11  as  one  (There  Is  no  l  -  m  *  0  spacecharge 

wave),  m  must  be  quite  large.  Consequently,  f±  always  Is  small  near  UQ 
*  0  for  spacecharge  modes,  and  resistive  wall  growth  Is  In  this  case  neg¬ 
ligible.  On  the  other  hand,  for  modest  positive  values  of  u0  Instability 
of  slow  spacecharge  waves  Is  described  by  Eq.  (14),  with  growth  rates  com¬ 
parable  to  those  of  slow  cyclotron  waves. 


Brelzman  and  Ryutov  have  shown  the  existence  of  another,  high  fre¬ 
quency  branch  of  the  resistive  wall  Instability  for  spacecharge 
modes. To  second  order  In  R/R0,  the  growth  rate  for  this  Insta¬ 
bility  Is  Identical  to  that  In  a  straight  tube.  Peak  growth  for  the  m  =  0 
mode  occurs  at  t  •  yR0/R,  with  y  the  normalized  electron  energy. 


(™) 


Note  that  this  expression  Is  Independent  of  the  magnetic  guide-field 
strength.  At  long  wavelengths  the  growth  rate  falls  off  as  ji1^,  until 
It  merges  with  the  lower  frequency  branch.  The  growth  rate  decreases  more 
gradually  at  shorter  wavelengths.  Interestingly,  the  Instability  persists 
for  m  #  0,  although  with  reduced  growth  rates.7 


>V»\V %•  •/  *.*  •:  os* 


Numerical  studies  with  the  laminar  beam  stability  code  GRADR^*^ 
confirm  the  absence  of  a  corresponding  high  frequency  cyclotron  wave  resis¬ 
tive  wall  Instability.7 

Me  Illustrate  the  relative  Importance  of  the  various  forms  of  the 
resistive  wall  Instability  for  the  typical  modified  betatron  parameters 
listed  In  Table  1.  A  stainless  steel  tube  wall,  with  normalized  conductiv¬ 
ity  o  *  5.24  •  10$  cur1.  Is  assumed.  (The  conductivity  of  copper  Is 
about  forty  times  larger.)  Inserting  these  values  Into  Eq  (II),  we  find 
r  *  4.0  •  10”5  carl  for  m  ■  1  cyclotron  waves  at  w0  ■  0.  The  corres¬ 
ponding  skin  depth  Is  £  *  0.07  cm;  the  cavity  wall  must  be  at  least  this 
thick  for  Eq.  (11)  to  be  valid.  The  growth  rate  value  just  given  assumes 
A  *  1.  The  Increase  of  t  near  the  end  of  the  acceleration  period  reduces  r 
by  less  than  25%.  Equation  (13)  Is  maximized  by  y  -  70,  where  the  t 
value  for  the  »0  -  0  resonance  drops  to  one.  Total  growth  is  rat  =  3.6 
in  this  case.  After  passing  through  the  w0  -  0  resonance,  the  i  *  m  *  1 
slow  cyclotron  mode  continues  to  grow,  at  the  rate  in  Eq.  (14).  Evaluating 
this  expression  for  the  typical  frequency  u0  -  Oz/y  yields  r  *  5.1  •  10"6  cm" 
Growth  occurs  over  about  1/3  the  acceleration  period,  or  rat  «  55  when  the 
<d0  «  0  contribution  Is  Included. 

As  already  discussed,  growth  of  slow  spacecharge  waves  near  u0  *  0 
Is  negligible.  At  higher  frequencies,  the  growth  rate  of  the  t  *  m  *  1 
spacecharge  wave  Is  comparable  to  that  of  the  corresponding  cyclotron 
mode.  Since  growth  occurs  throughout  the  acceleration  period,  rat  *  150. 

At  sti'1  higher  frequencies  m  *  0  spacecharge  waves  grow  at  a  maximum  rate 
of  r  *  1.3  •  10-6,  or  about  25%  of  the  Intermediate  frequency  growth 
rate.  It  Is  Important  to  bear  In  mind  that  these  growth  estimates  probably 
are  uncertain  by  a  factor  of  two,  due  to  approximations  made. 
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TABLE  1 


Typical  Modified  Betatron  Parameters  Used  In  Evaluating  Resistive  Wall 

Instability  Growth  Rates. 


Toroidal  Magnetic  Field 
Vertical  Magnetic  Field  (Initial) 
Toroid  Major  Radius 

Toroid  Minor  Radius 

Beam  Radius 

Beam  Current 

Beam  Energy 


2.5  kg 
115  g 
100  cm 

10  cm 

1  cm 

0.59 

7-100 


Acceleration  Time 


T 


3.107  cm 


The  preceding  analysis  Ignores  thermal  spread  In  the  beam  electron 
velocities.  As  a  result  of  the  toroidal  geometry,  a  spread  In  electron 
toroidal  velocity  can  smear  out  the  transverse  resonances  associated  with 
the  resistive  Instabilities.  In  a  companion  paper  we  find  that  the 
required  spreads  In  toroidal  and  transverse  velocities  are  achievable  from 
a  few  percent  spread  In  Initial  electron  kinetic  energy.14  A  similar 
conclusion  was  drawn  In  Ref.  9,  based  on  a  Lorentzlan  distribution  of  ener 
gles.  Of  course,  this  damping  of  spacecharge  and  cyclotron  wave  Instabil¬ 
ities  by  a  spread  In  electron  velocities  Is  Ineffective  for  low  frequency 
modes.  As  we  have  already  seen,  however,  these  latter  modes  are  not  a  ser 
lous  problem. 

In  summary,  we  find  that  the  m  ■  1,  low  i  spacecharge  and  cyclotron 
wave  Instabilities  are  the  most  dangerous  of  the  various  resistive  wall 
phenomena  Identified  for  modified  betatrons.  Amplification  factors  of 
e150  -  1065  and  e55  -  1024,  respectively,  are  predicted  for 
the  parameters  In  Table  1.  Successful  operation  of  a  modified  betatron 
requires  cutting  the  growth  exponent  rat  to  about  unity.  (Computer  Simula 
tlons  suggest  that  Initial  perturbations  may  be  quite  large.15) 

Employing  more  highly  conducting  cavity  walls  provides  a  factor  of  six. 
Reducing  the  acceleration  time  by  an  order  of  magnitude  would  then  effect¬ 
ively  eliminate  these  Instabilities.  Alternatively,  an  Initial  spread  of 
few  percent  In  the  electron  energy  should  be  sufficient  to  damp  out  the 
modes. 

We  ere  Indebted  to  M.  Campbell,  M.  Jones,  D.  Sullivan,  and  particu¬ 
larly  P.  Sprangle  for  helpful  discussions.  This  research  Is  supported  by 
Office  of  Naval  Research  contract  N00014-81-C-0647,  monitored  by  C. 
Roberson. 
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I.  INTRODUCTION 


High  current  racetrack  beam  Induction  accelerators  and  modified  beta¬ 
trons  are  a  subject  of  Increasing  Interest  as  sources  of  high  power  elec¬ 
tron  beams  for  free  electron  lasers,  flash  radiography,  and  other  applica¬ 
tions.  The  racetrack  Induction  accelerator  geometry  is  illustrated  in  Fig¬ 
ure  1.  The  beam  Is  Injected  from  a  conventional  pulsed  diode  beam  genera¬ 
tor  Into  the  drlfttube.  Is  progressively  accelerated  as  it  repetitively 
passes  one  or  more  Induction  modules,  and  then  is  extracted  from  the  accel¬ 
erator  for  Its  Intended  use.  Beam  extraction  may  even  be  unnecessary  for 
microwave  applications,  because  a  slow-wave  or  rippled-magnetic-field  cav¬ 
ity  can  be  Inserted  In  a  straight  section  of  the  drlfttube.1 

Host  beam  stability  studies  for  high  current  recirculating  devices 
have  dealt  with  negative  mass  and  resistive  wall  Instabilities.2"5  How¬ 
ever,  experience  with  linear  induction  accelerators  suggests  that  beam 
breakup  and  Image  displacement  Instabilities  due  to  beam  interaction  with 
the  Induction  modules  and  other  discontinuities  In  the  drifttube  may  be 
significant.6*7  The  beam  breakup  Instability  arises  from  a  resonant 
coupling  between  beam  transverse  oscillations  and  m=l  electromagnetic  cav¬ 
ity  modes  localized  to  the  acceleration  gaps,8  while  the  Image  displacement 
instability  Is  caused  by  interrupting  the  m*l  beam  Image  current  in  the 
drlfttube  wall.9  Clearly,  the  two  are  Interrelated;  in  some  contexts  the 
image  displacement  instability  can  be  viewed  as  the  low  frequency  limit  of 
the  beam  breakup  instability. 

This  paper  extends  Instability  results  developed  for  high  current  lin¬ 
ear  induction  accelerators  to  cyclic  devices.  The  primary  differences  are 
(1)  that  the  beam  passes  the  same  few  gaps  again  and  again,  and  (2)  that 
the  curved  beam  trajectory  couples  longitudinal  to  transverse  motion,  per¬ 
haps  leading  to  a  hybridization  of  the  negative  mass  and  beam  breakup 
modes.  Experience  with  Interaction  among  Instabilities  In  other  situations 
suggests  that  the  latter  Item  Is  the  less  Important  of  the  two.10 


Therefore,  we  consider  only  the  former  aspect,  periodicity  of  the  inter¬ 
action,  and  omit  curvature  effects.  As  In  earlier  calculations,  we  assume 
that  the  beam  propagates  at  the  speed  of  light,  that  lateral  displacements 
of  the  beam  are  small  and  rigid,  and  that  the  gaps  are  narrow  compared  to 
the  beam  oscillation  wavelength.  The  consequences  of  relaxing  these 
approximations  for  the  Image  displacement  Instability  under  different  cir¬ 
cumstances  Is  discussed  elsewhere.9  We  expect  the  three  approximations  to 
be  adequate  for  present  purposes. 

A  general  dispersion  relation  for  a  multiple  gap  racetrack  Induction 
accelerator  Is  derived  In  Section  II.  Because  the  resulting  expression,  a 
determinant.  Is  cumbersome,  we  then  specialize  to  the  case  of  Identical, 
uniformly  spaced  gaps,  obtaining  approximate  analytical  growth  rate  form¬ 
ulas  for  the  beam  breakup  and  Image  displacement  Instabilities.  Section 
III  illustrates  Instability  growth  for  a  1  kA  electron  beam  accelerated  to 
40  MeV  through  fifty  cycles  of  a  four  gap  device.  (These  parameters  are 
based  on  an  Induction  llnac  developed  by  the  National  Bureau  of  Standards 
and  now  being  modified  at  the  Naval  Research  Laboratory.11)  First,  growth 
rates  are  estimated  based  on  the  formulas  of  the  preceding  Section.  Then, 
more  precise  numerical  values  are  obtained  with  the  beam  transport  code 
BALTIC,  including  the  effects  of  acceleration  and  initial  transients.  Ou? 
findings,  summarized  In  Section  IV,  are  encouraging:  The  image  displace¬ 
ment  Instability  Is  negligible,  and  the  beam  breakup  Instability  is  manage¬ 
able  provided  that  the  mode  quality  factor  Q  Is  kept  low.  A  comparison 
with  the  negative  mass  and  resistive  wall  Instability  growth  rates  also  is 
provided. 


II.  DERIVATION  OF  GROWTH  RATE  EXPRESSIONS 

In  the  long  wavelength,  paraxial  approximation  the  linearized  equation 
of  transverse  motion  for  the  beam  centroid  is  given  by9*10 


Y  »,^ 
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Here,  e  «  x  +  1y  Is  the  (complex)  transverse  displacement  of  the  beam 
centroid,  p  is  the  beam  density,  y  Is  the  beam  relativistic  energy,  a  is 
the  beam  radius,  R  is  the  drifttube  radius,  and  «c  is  the  cyclotron 
frequency  for  the  magnetic  guide  field.  Time  is  normalized  to  the  beam 
axial  velocity,  which  Is  in  turn  assumed  to  be  equal  to  the  speed  of 
light.  The  last  term  on  the  left  side  of  Eq.  (1)  represents  the  combined 
transverse  forces  on  the  beam  from  Its  Image  charge  and  current  In  the 
smooth  drifttube  wall.  Discontinuities  In  the  wall  are  treated  as 
impulsive  forces  appearing  on  the  right  side  of  Eq.  (1). 


We  begin  the  derivation  of  the  dispersion  matrix  by  assuming  the  beam 
parameters  to  be  constant  in  time  and  Fourier  transforming  Eq.  (1). 
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As  mentioned  In  the  Introduction,  we  shall  take  account  of  time-varying  y 
In  Sec.  Ill  by  direct  numerical  Integration  of  Eq.  (1). 

Next,  Eq.  (2)  Is  solved  for  5,  treating  the  right  side  of  the  equation 
as  a  source  term.  This  can  be  done  using  a  Green's  function  or,  equiva¬ 
lently,  by  twice  Integrating  the  equation  and  Imposing  periodicity  at  z  * 
L,  the  accelerator  path  length.  The  result  is 
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C(Zj )  exp[^+(si9n[2j  -  z]L/2  -  (z..  -  z))] 
C(Zj)  exp[ik_(s1gn[Zj  -  z]L/2  -  (z^  -  z))] 


where  N  Is  the  number  of  acceleration  gaps  and  S1gn(x)  Is  the  sign  of  x. 
The  wavenumbers  It*  characterize  the  beam  cyclotron  and  spacecharge  waves, 
respectively,  as  determined  by  the  left  side  of  Eq.  (2),  i.e.,  between 
gaps. 
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For  typical  Induction  linac  parameters,  k+  -  u  +  u»c/Y  and  k_  -  u. 
Although  It  may  t  be  Immediately  apparent,  e(z)  as  expressed  In  Eq.  (3) 
Is  everywhere  continuous;  Its  first  derivative  is  discontinuous  at  Zj. 


Evaluating  Eq.  (3)  at  each  of  the  gaps  yields  the  system  of  equations 
v(»<+  -  O  C  (zA)  « 

(5) 

(2s1nk+L/2)-1  J  Fj  c(Zj)  exp[1k+(L/2  -  Mod[Zj  -  z£,L])] 

-  (2s1nk_L/2)-1  J  Fj  c(2j)  exp[1k_(l/2  -  Mod[Zj  -  z£,L])] 

J*1 

Here,  Hod  (x,L)  Is  x  modulo  L,  defined  to  lie  between  0  and  L.  The  deter¬ 
minant  of  the  coefficient  matrix  of  Eq.  (5)  Is  the  desired  dispersion  rela¬ 
tion.  In  general.  It  must  be  solved  numerically  to  obtain  «.  Additional 
analytical  progress  can,  however,  be  made  when  the  gaps  are  spaced  uni¬ 
formly  and  all  have  the  same  response  function  F. 
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Y(k+  -  kj  C  (tL/N) 


(6) 


COS[ (k+  -  k_)L/2N] 


»  "  « 


>  V 


For  F/**c  small  those  equation  are  approximately  satisfied  by  k+  -  kn  or 
k.  -  kn,  with  kn  a  2m/L.  Correspondingly. 

«  -  *  s  k  -  («  ±  [»,2  -  2Pa2/vR2]1/2)/2v  •  U°) 

n  n  v  c  c 

Instability  may  occur  when  u  matches  a  resonant  frequency  of  F  or  when  the 
right  side  of  Eq.  (9)  exceeds  unity  away  from  a  resonance.  These  two  pos¬ 
sibilities  are  realized  In  the  beam  breakup  and  image  displacement  Insta¬ 
bilities,  respectively. 


The  beam  breakup  Instability  Is  caused  by  TMln0  oscillatory  elec¬ 
tromagnetic  fields  localized  to  the  acceleration  gaps  and  adjacent  drift- 
tube  regions.  These  cavity  modes  are  conveniently  represented  by  damped 
harmonic  oscillator  equations  for  their  normalized  vector  potentials. 


<*>- 


TT 
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e(z1 ) 


(ii) 


Q  Is  the  cavity  mode  quality  factor,  while  w0  <*  its  frequency.  Coupling 
between  the  beam  and  a  cavity  mode  Is  given  by  Its  transverse  Impedance 
Zj/Q.  The  transverse  Impedance  Is  basically  a  geometrical  factor,  which 
tends  to  scale  linearly  with  the  gap  width. 


Fourier  transforming  Eq.  (11)  In  time,  we  obtain  for  F  -  A/£ 


F 


♦  lUWg/Q  “  Up 
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(12) 
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The  beam  breakup  Instability  growth  rate  now  can  be  approximated  In  the 
usual  manner  by  expanding  Eq.  (8)  and  (12)  about  «n,  defined  In  Eq. 
(10).  The  resulting  quadratic  In  u  -  un  shows  Instability  only  for  kn 
■  k-»  with  growth  rate 


Note  that  this  derivation  falls  for  s1n[(k+  -  k_)L/2]  too  small;  see  Eq. 
(9).  More  precisely,  the  Instability  growth  rate  Is  reduced  whenever 
| Nkn  -  wc/v|  Is  less  than  the  absolute  value  of  the  right  side  of  Eq. 

(13)  and  finally  vanishes  when  Nkn  -  wc/y. 

The  Image  displacement  Instability  arises  due  to  Interruption  of  the 
beam  Image  current  at  a  discontinuity  In  the  drlfttube  wall,  such  as  occurs 
at  an  acceleration  gap.  The  analysis  In  Ref.  9  and  12  gives  a  force  coef¬ 
ficient 

F  -  K  eT-  (H) 

The  effective  gap  width  *  Is  given  by  the  physical  width  when  the  latter  is 
small  compared  to  the  drlfttube  radius  R.  In  the  opposite  limit,  i  <  R. 
As  already  noted,  growth  occurs  whenever  the  right  side  of  Eq.  (9)  exceeds 
unity.  For  F/a»c  small  this  happens  only  In  a  narrow  parameter  region 
centered  on  tan[(k+  -  k„)l/2]  »  F/y(k+  -  k_),  for  which  the  growth  rate  Is 


r  •  MF/tucL 


(15) 


In  the  next  section  we  evaluate  these  two  growth  rates  for  possible 
recirculating  Induction  accelerator  parameters  and  then  generalize  our  num¬ 
erical  results  to  Include  Irregular  gap  spacing  and  beam  acceleration  by 
computer  solution  of  Eq.  (1). 

III.  NUMERICAL  EVALUATION  OF  GROWTH  RATES 

The  four  gap  linear  Induction  accelerator  at  the  Naval  Research  Labor¬ 
atory  Is  an  attractive  candidate  for  a  recirculating  device.  The  accelera¬ 
tion  gaps  are  In  pairs  with  the  members  separated  by  about  30  cm.  The 
distance  between  the  two  pairs  Is  arbitrary  within  reasonable  limits.  Here 
we  take  the  distance  to  be  200  cm,  giving  a  total  round-trip  path  length  of 
L  ■  460  cm.  (The  specific  value  chosen  does  not  strongly  Influence  our 
conclusions.)  The  principle  m*l  gap  normal  mode  Is  at  880  MHz  with  a 
transverse  Impedance  of  15  ohms  and  a  Q  of  60.  Each  gap  provides  an  accel¬ 
eration  of  0.2  MeV,  resulting  In  an  energy  gain  for  the  electron  beam  of  40 
MeV  after  50  passes.  A  beam  current  of  I  kA  and  guide  field  of  2  kg  are 
assumed.  These  parameters  are  summarized  In  Table  1.  Note  that  y  In  the 
table  has  been  corrected  for  the  spacecharge  depression  In  the 

drlfttube  and  also  that  Zj/Q  has  been  expressed  In  dimensionless  form  by 
dividing  the  transverse  impedance  by  30  ohms. 

Although  the  four  gaps  are  not  spaced  uniformly,  a  reasonable  upper 
bound  on  beam  breakup  Instability  growth  can  nonetheless  be  obtained  from 
Eq.  (13).  We  find  r  »  12.9*  10-**  cm-1,  corresponding  to  about  30 
e-fol dings  during  the  course  of  acceleration.  (We  assume  wo  ■  wn.) 
The  most  obvious  way  of  cutting  this  amplification  to  an  acceptable  level 
Is  hy  reducing  Q.  Recent  work  at  Lawrence  Livermore  National  Laboratory  on 
short  pulse  Induction  accelerators  suggests  that  a  quality  factor  as  low  as 


TABLE  1.  NOMINAL  RACETRACK  INDUCTION  ACCELERATOR  PARAMETERS 


Path  Lengths 
Drlfttube  Radius 
Beam  Radius 
Guide  Field 
Beam  Current 
Beam  Energy 
Number  of  Revolutions 
Number  of  Gaps 
Acceleration  per  Gap 
Gap  Resonant  Frequency 
Mode  Quality  Factor 
Gap  Transverse  Impedance 
Gap  Width 


L  *  460  cm 
R  ■  7  cm 
a  «  1  cm 

Bz  ■  2  kg  (wc  *  1.173  cm' 

1  «  1  kA  ( v  *  0.0588) 

U  -  0.4  -  40  MeV  (y  -  1.5  -  80) 
50 

N  =  4 

AU  *  0.2  MeV  (Ay  -  0.4) 

880  MHz 
Q  «  60 
15  ohms 
i  ■  5  cm 


(Zi/Q  -  0.5) 


six  Is  achievable.13  For  Q-6,  r  •  2.4»10-4  cm-1,  or  5.5  e-foldings.  As 
many  as  eight  e-fol dings  may  be  tolerable  provided  the  gaps  are  not  excited 
appreciably  a:  ‘heir  resonant  frequency  prior  to  beam  Injection. 

More  precise  growth  rates  can  be  obtained  by  numerically  solving  Eq. 
(5)  plus  Eq.  (12),  and  this  Is  not  difficult  to  do.  Instead,  we  choose  to 
Integrate  Eq.  (1)  plus  Eq.  (11)  directly  using  the  computer  program 
BALTIC.  The  code  was  exercised  extensively  In  support  of  the  RADLAC  radial 
pulsellne  accelerator  program  and  so  Is  well  tested.10  BALTIC  has  the 
advantage  over  a  dispersion  relation  solver  that  It  takes  account  of  beam 
acceleration  and  of  transients.  It  Is,  of  course,  much  slower. 

Table  2  summarizes  thirteen  runs  of  the  BALTIC  code.  The  quantities 
N,  0,  u>0,  and  y  were  varied,  but  with  NZ^/Q  held  fixed.  Some  calculations 
involved  an  accelerating  beam.  Cases  1-4  illustrate  the  effect  of  changing 
u0  -  «n*  Because  kn  «  0.014-n  cm-1,  the  pattern  of  growth  rate 
variation  with  Uo  In  these  four  cases  repeats  Itself  with  the  same  0.014 
cm-1  periodicity  In  «0.  It  so  happens  that  s1n[(k+  -  k_)L/2]  is  small 
for  y  ■  80,  the  value  chosen  In  cases  1-4,  which  reduces  the  growth  rate 
some.  Setting  y  *  60  avoids  this  situation,  increasing  r  by  25%,  as 

Illustrated  In  cases  5  and  6.  The  growth  rate  change  caused  by  dropping  Q 

from  60  to  10,  as  In  case  7,  Is  consistent  with  Eq.  (13). 

The  seven  cases  just  discussed  were  for  a  single  gap.  Cases  8  and  9 

treat  two  gaps  evenly  spaced.  Growth  rates  are  comparable  to  the  single 
gap  runs,  as  expected.  The  enhancement  of  r  for  case  9  relative  to  case  8 
again  arises  from  adjusting  s1n[(k+  -  k.)L/2], 

Finally,  cases  10-13  treat  four  gaps  spaced  as  described  at  the  begin¬ 
ning  of  this  Section.  Also,  the  beam  Is  accelerating  In  these  last  four 
runs.  It  Is  clear  from  the  resulting  growth  rates  that  the  gradual  accel¬ 
eration  which  occurred  has  little  effect  on  stability  except  to  average 
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TABLE  2.  SUMMARY  OF  BEAM  BREAKUP  INSTABILITY  CALCULATIONS  WITH  BALTIC 


Case 

N 

_Q 

“0, 

_  (cm-1 ) 

Y 

r 

(10-4  cm 

1 

1 

60 

0.18 

80 

8.3 

2 

1 

60 

0.1732 

80 

3.3 

3 

1 

60 

0.17757 

80 

8.9 

4 

1 

60 

0.1787 

80 

8.8 

5 

1 

60 

0.18 

60 

10.0 

6 

1 

60 

0.17757 

60 

12.3 

7 

1 

10 

0.18 

80 

3.4 

8 

2 

60 

0.18 

80 

9.6 

9 

2 

60 

0.18 

1.5 

12.0 

10 

4 

60 

0.18 

1.5-80 

9.8 

11 

4 

10 

0.18 

1.5-80 

2.6 

12 

4 

6 

0.18 

1.5-80 

1.2 

13 

4 

6 

0.17757 

1.5-80 

1.6 

over  n>0  -  un  and  s1n[(k+  -  k_)L/2].  In  case  13,  w0  was  decreased 
slightly  to  maximize  growth  at  large  -y.  Indeed,  r  Increased  by  one-third 
relative  to  case  12. 

Figure  2  depicts  the  maximum  transverse  displacement  of  the  beam  as  a 
function  of  time  for  case  13.  It  and  the  other  runs  were  Initiated  with  a 
uniform  transverse  offset  of  unit  magnitude  for  the  beam.  The  cavity  modes 
were  Initially  excited  at  an  amplitude  consistent  with  Injection  of  the 
offset  beam  with  zero  risetime,  a  worst  case  assumption.  The  displacement 
Is  seen  to  grow  by  a  factor  of  about  3.6. 

The  Image  displacement  Instability  turns  out  to  be  Insignificant. 
From  Eq.  (15)  its  peak  growth  rate  Is  0.9»10-4  cm-1.  Since  this  occurs 
only  over  a  narrow  range  of  parameters,  however,  we  should  expect  a  much 
slow,  r  average  growth  during  acceleration.  Indeed,  a  four  gap  BALTIC  run 
gave  no  coherent  growth  whatsoever,  and  transient  fluctuations  amounted  to 
less  than  15%  of  the  Initial  displacement. 

IV.  CONCLUSIONS 

In  this  report  we  have  developed  a  simple  theory  of  bean  breakup  and 
Image  displacement  Instabilities  in  cyclic  Induction  accelerators  and 
applied  It  to  obtain  growth  rate  estimates  for  a  possible  device.  The 
theoretical  model  takes  recirculation  Into  account  primarily  by  enforcing 
periodicity  on  the  unstable  transverse  models.  Drifttube  curvature  effects 
were  not  considered.  Me  find  for  the  1  kA,  40  MeV  accelerator  that  neither 
the  beam  breakup  nor  the  Image  displacement  Instability  should  be  serious. 
Specifically,  with  and  parameters  otherwise  as  In  Table  1  the  bean 
breakup  growth  rate  Is  r  «  1.6* 10-4  cm-1;  the  Image  displacement  growth  is 
negligible. 


13 


S'a'Mv! 


For  comparison,  growth  rates  of  the  negative  mass  and  resistive  wall 
(for  stainless  steel)  Instabilities  for  comparable  parameters  are  of  order 
5»10-4  cm*1  and  0.2«10-4  cm-1,  respectively.3-5  Thus,  the  negative  mass 
Instability  Is  somewhat  more  serious  than  the  gap-induced  instabilities  in 
this  Instance. 

The  accelerator  path  length  L  *  460  cm  chosen  for  the  numerical  calcu¬ 
lations  Is  perhaps  small  for  practical  devices.  Increasing  L  by  lengthen¬ 
ing  the  straight  sections  of  the  racetrack  geometry  would  have  little 
1mpart  on  any  but  the  resistive  wall  instability  integrated  growth:  Growth 
rates  of  the  gap-induced  instabilities  decrease  roughly  as  L-1,  so  that 
growth  per  pass  remains  fixed,  while  the  negative  mass  instability  grows 
only  In  the  curved  sections  of  the  drift  tube.  In  fact,  the  negative  mass 
mode  may  decay  in  straight  sections.  Growth  per  cycle  of  the  resistive 
wall  Instability,  of  course,  increases  linearly  with  L.  Its  growth  rate  is 
so  slow,  however,  that  the  path  length  could  be  increased  by  an  order  of 
magnitude  or  more  before  the  resistive  wall  instability  became  competitive 
with  the  other  modes. 

The  Important  stabilizing  influences  of  electron  energy  spread  and 
transverse  oscillations  have  been  omitted  from  our  calculations.  Recent 
estimates  for  betatrons  with  comparable  parameters  suggest  that  a  10% 
spread  in  the  injected  beam  electron  energies  is  sufficient  to  suppress  the 
negative  mass  Instability  provided  that  the  mean  injection  energy  is 
greater  than  about  3  MeV.2-5  Additionally,  research  in  support  of  the  ATA 
accelerator  program  predicts  that  employing  occasional  nonlinear  focusing 
elements  to  spread  the  electron  transverse  oscillation  frequencies  can 
limit  beam  breakup  instability  growth.14  He  remark  also  that  varying  the 
resonant  frequencies  h>0  From  gap  to  gap  can  cut  the  beam  breakup  growth 
rate  by  N-1 ,  provided  the  relative  shifts  are  greater  than  Q-1.15 
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Extensions  to  the  present  research  naturally  fall  Into  two  categories, 
adding  drlfttube  curvature  and  adding  beam  temperature.  The  former  is  rel¬ 
atively  easily  accomplished  by  inserting  the  transport  equations  of  Ref.  4 
or  16  Into  BALTIC.  This  also  would  allow  us  to  determine  the  Influence  of 
straight  drlfttube  sections  on  the  negative  mass  instability.  Treating 
beam  temperature  beyond  what  has  already  been  done  probably  requires  numer¬ 
ical  simulation. 
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SUMMARY 


Beam  breakup  and  negative  mass  Instability  growth  rates  for  a  1  kA,  40 
MeV  electron  beam  racetrack  Induction  accelerator  are  computed.  The  device 
9  1$  taken  to  have  four  acceleration  gaps,  each  with  0.2  MeV  applied  voltage 

and  15  ohm  transverse  impedence;  the  guide  field  Is  2  kg.  We  find  that  the 
total  amplification  of  the  beam  breakup  mode  Is  limited  to  five  e-foldings 
provided  that  the  cavity  mode  quality  factor  Q  Is  6.  Thus,  the  negative 
0  mass  Instability,  which  grows  several  times  faster,  Is  the  dominant  con¬ 
sideration.  However,  we  also  find  that  the  energy  range  over  which  the 
negative  mass  Instability  occurs  can  be  narrowed  substantially  by  reducing 
the  guide  field  strength  after  the  beam  has  been  accelerated  to  about  12 
I  MeV.  This  approach,  coupled  with  beam  thermal  effects,  not  considered 

here,  probably  is  sufficient  to  limit  negative  mass  growth  to  acceptable 
levels  in  the  racetrack  accelerator. 

INTRODUCTION 

High  current  racetrack  induction  accelerators  and  modified  betatrons 
are  a  subject  of  increasing  Interest  as  sources  of  high  power  electron 
beams  for  free  electron  lasers,  flash  radiography,  and  other  applications. 
The  racetrack  Induction  accelerator  geometry  is  illustrated  schematically 
in  Figure  1.  The  beam  is  Injected  from  a  conventional  pulsed  diode  bean 
generator  Into  the  drlfttube,  is  progressively  accelerated  as  it  repeti¬ 
tively  passes  one  or  more  induction  modules,  and  then  Is  extracted  from  the 
accelerator  for  Its  Intended  use.  Extraction  may  even  be  unnecessary  for 
microwave  applications,  because  a  slow-wave  or  rippled-magnetic-field  cav¬ 
ity  can  be  Inserted  In  a  straight  section  of  the  drifttube.1 

Most  beam  stability  studies  for  high  current  recirculating  devices 
have  dealt  with  the  negative  mass  and  resistive  wall  instabilities.2'5 
However,  experience  with  linear  Induction  accelerators  suggests  that  beam 
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breakup  due  to  Interaction  with  the  Induction  modules  and  other  discontinu¬ 
ities  In  the  drlfttube  may  be  significant.6*7  The  beam  breakup  instabil¬ 
ity  arises  from  a  resonant  coupling  between  beam  transverse  oscillations 
and  m»l  electromagnetic  cavity  modes  localized  to  the  acceleration  gaps, 
resulting  In  large  lateral  displacements  of  the  beam.6*9  In  this  pape- 
we  present  a  linear  dispersion  relation  describing  both  beam  breakup  and 
negative  mass  Instabilities,  Including  their  possible  interaction,  and 
evaluate  it  for  parameters  of  the  proposed  racetrack  induction  accelerator 
designed  by  the  Naval  Research  Laboratory.2 

The  NRL  device  is  based  on  the  four  module  linear  induction  acceler¬ 
ator  developed  by  the  National  Bureau  of  Standards.10  It  is  expected  to 
accelerate  a  1  kA  electron  beam  from  1  to  40  MeV  in  fifty  cycles.  The  beam 
and  drifttube  radii  are  1  and  7  cm,  respectively.  The  principle  m=l  reson¬ 
ance  of  the  gaps  has  a  frequency  of  880  MHz,  an  impedence  of  15  ohms,  and 
a  quality  factor  (Q)  of  60.  Experience  with  the  ETA  linear  Induction 
accelerator  at  Lawrence  Livermore  National  Laboratory  indicates  that  Q  can 
be  greatly  reduced,  however,  and  we  shall  take  0*6  in  our  numerical  work.11 
The  NRL  design  includes  a  2  kg  axial  magnetic  field  to  maintain  the  beam 
equilibrium  and  Improve  beam  stability  at  low  energies.  Reducing  or  elim¬ 
inating  this  guide  field  at  higher  energies  is  nonetheless  an  interesting 
possibility.  The  stability  analyses  below  consider  both  options. 

DISPERSION  RELATION 

For  simplicity  we  represent  the  racetrack  accelerator  as  a  torus  with 
a  single  gap.  These  two  approximations  are  conservative  in  that  omitting 
the  straight  sections  of  the  racetrack  and  lumping  the  several  gaps  into 
one  overestimate  negative  mass  and  beam  breakup  growth,  respectively.  The 
desired  dispersion  relation  is 


(fl2-«r2  ♦  FrM  +  we2Y2x)  (»2-tt>22  +  F2/yL) 


-  (Bft/Y)2  fl2  -  0 


x  ■  v  (uv0  -  ^)/[R2  +  v  (w2 
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4)3 

IT 


v  *  (1  +  2m  a/r^)  v/y' 


Here,  n  ■  w-luig  is  the  Doppler-shifted  wave  frequency,  i  is  the  toroidal 
mode  number,  and  u>8  *  vg/R  is  the  toroidal  rotation  frequency  of  the 
beam.  (The  poloidal  rotation  frequency  is  assumed  negligible.)  The  radial 
and  vertical  betatron  frequencies  are 

2m  2,  2»  2 

«r  *  (1  -  n  -  ns  rfa  /a  )  we 


("  -  "s  rbZ''a2)  V 


"s  5  "t./(2“e2',3) 


with  n  the  betatron  index,  n^  the  beam  density,  v  =  %rt)2/2 
Budker's  parameter,  and  y  the  beam  energy.  The  drifttube  major  radius  is 
R,  the  drifttube  minor  radius  is  a,  and  the  beam  minor  radius  is  r^. 
L«2nR.  The  toroidal  guide  field  strength  is  B9;  the  betatron  field 
strength  enters  as  Bz  ■  -wgy. 
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In  a  high  current  betatron  «r2  and  w22  can  be  of  either  sign.  The 
beam  Is  unstable,  however,  whenever 

»g2  *  »22  wr2/(80/Y)2  (5) 

Is  negative.  To  avoid  this  situation,  as  well  as  for  simplicity,  we  take 
n*l/2.  The  energy  at  which  og2  «  0  typically  is  labeled  the  transition 
energy, 

ytr  -  (4vR2/a2)1/3  (6) 


The  gap  response  function  F  is  defined  as9 


where  o>0  is  the  resonant  frequency,  Z^/Q  is  the  transverse  impedence,  and 
Q  is  the  quality  factor.  Setting  F=0  in  (1)  recovers  the  high  current  bean 
negative  mass  dispersion  relation.  The  negative  mass  instability  occurs 
for  all  t  over  a  broad  range  of  energies  when  y  >  ytp.  For  low  t  only, 
one  or  two  instability  bands  (often  overlapping)  also  may  exist  when  y  < 
Ytr*  Three  the  s1x  beam  modes  (m»0  spacecharge,  m=l  spacecharye, 
and  m*l  cyclotron;  m  is  the  poloidal  mode  number)  have  negative  energy  and 
so  can  couple  unstably  to  the  gap  fields.  Note  that  coupling  in  the  m*0 
spacecharge  mode  occurs  only  due  to  toroidal  curvature.  Choosing  R=70  cm, 
we  find  maximum  coupling  at  i»13. 
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LARGE  B  ANALYSIS 


For  the  parameters  considered  here  and  toroidal  mode  numbers  In  the 
vicinity  of  13,  the  negative  mass  Instability  exists  only  beyond  the  trans¬ 
ition  energy  -ytr  ■  2.9.  Just  above  the  transition  energy  the  insta¬ 
bility  Is  due  solely  to  the  Interaction  between  the  positive  and  negative 
energy  m-1  spacecharge  modes,  while  at  still  higher  energies  the  m=0  space- 
charge  modes  also  are  Involved.  This  change  is  readily  visible  in  the  neg¬ 
ative  mass  Instability  growth  rate,  the  dashed  curve  in  Figure  2.  Which 
portion  grows  faster  depends  on  circumstances. 

Although  the  peak  growth  rate  at  lower  energy  is  not  readily  deter¬ 
mined  analytically,  the  higher  energy  peak  is  easily  shown  to  be 


/3  ro  2  2-.1/3  /o\ 

r  •  -jn  [2  l  u>0  u>g  v  y  J 

Instability  ceases  for 

y  >  [6/3  *RB0  V  (1  ♦  2  in  a/rb)31/2  (9) 


here  about  62. 


In  the  absense  of  curvature,  the  beam  breakup  growth  rate  also  is 
easily  estimated.  For  Q  not  too  large,9 


v  Z^/Q 

r  -  “oQ  -B^r 


(10) 


Both  m*l  negative  energy  modes  grow  at  this  rate  when  their  frequencies 
roughly  match  u0. 
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The  solid  curves  In  Figure  2  show  growth  rates  of  the  negative  mass 
and  beam  breakup  Instabilities  combined.  The  negative  mass  results  are 
seen  to  be  only  weakly  affected  by  the  gap  resonance.  The  m=l  cyclotron 
and  hybrid  m*0/l  spacecharge  modes  have  become  unstable,  however,  with  a 
growth  rate  agreeing  with  (10)  to  within  a  factor  of  1.5.  These  findings 
are  insensitive  to  small  changes  in  the  resonant  frequency. 

SHALL  B  ANALYSIS 


Although  modified  betatron  and  racetrack  induction  accel  rat or  studies 
usually  assume  a  large  toroidal  guide  field,  large  Be  is  in  fact  needed 
to  provide  a  beam  equilibrium  only  for  y  small.  Reducing  or  perhaps  elim¬ 
inating  Be  after  the  beam  has  been  accelerated  sufficiently  has  certain 
advantages  for  stability,  as  we  see  below. 


For  Bq=0,  the  negative  mass  growth  rate  is  approximately 

_  VT  ra  3  *  2-.1/3 

r  -  “2  L£  O)0  V  y  J 


(11) 


This  expression  exceeds  (8)  whenever  Be/y  >  u»Q.  However,  the  comes 
ponding  high  energy  cutoff. 


y  >  3/E  tv  (1+  2  tn  a/r&)  (12) 

here  27.5,  typically  is  much  lower  than  (9).  See  the  dashed  curve  in 
Figure  3. 


The  beam  breakup  instability  maximum  growth  rate  is  again  readily 
estimated,  this  time  giving 

v  Zj/Q 

r  *  “o  Q  'ZufC  y-  (1J) 
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Equation  (13)  exceeds  (10)  for  Bg/y  >  2 The  solid  curves  in  Figure 
3  show  the  effects  of  F*0.  As  in  Figure  2,  the  negative  mass  instability 
is  only  slightly  modified  by  the  gap;  the  beam  breakup  instability  is  des¬ 
cribed  reasonably  well  by  (13). 

A  comparison  of  the  two  figures  suggests  that  some  reduction  in  total 
Instability  growth  during  acceleration  can  be  achieved  by  rapidly  decreas¬ 
ing  the  guide  field  as  the  beam  energy  exceeds  about  12  MeV. 


RECOVERING  "CONVENTIONAL"  NEGATIVE  MASS  INSTABILITY  BEHAVIOR  FROM  THE 
HIGH-CURRENT-BETATRON  DISPERSION  RELATION 


The  preceding  small  Be  analysis  predicts  a  high  energy  cutoff  for 
the  negative  mass  instability  in  a  high  current  betatron,  while  standard 
derivations  of  negative  mass  growth,  performed  for  a  low  current  betatron, 
lead  to  no  such  cutoff. 12»13The  source  of  this  apparent  discrepancy, 
namely  the  failure  of  the  usual  approximation 

1  -  uj2R2/*2  -  y”2  (14) 

when  vy4  becomes  large,  has  been  noted  previously.4  To  emphasize  this 
point,  we  here  recover  the  accepted  low  current  growth  rate  and  show  that 
it  too  exhibits  a  high  energy  cutoff,  where  (14)  breaks  down.  Interest¬ 
ingly,  a  comparison  of  the  new  cutoff  with  (12),  the  high  current  limit, 
indicates  that  practical  betatron  parameters  exist  for  which  the  negative 
mass  instability  does  not  occur  at  all  at  moderate  toroidal  mode  numbers. 


We  begin  by  setting  Bg  (and  F)  to  zero  in  (1), 


n2  2,2  2  K  n 

0  -  u>r  +  u>0  y  x 


Next,  we  drop  n2  as  "'ll  compared  to  uy2, valid  for  m*0  spacecharge  waves 
in  low  current  beams,  and  rearrange  terms  in  the  expression  for  x* 
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This  derivation  shows  clearly  the  the  negative  mass  instability  at  low  cur¬ 
rents  Is  associated  only  with  the  m«0  spacecharge  waves.  At  higher  cur¬ 
rents*  for  which  ft  and  up  become  comparable,  the  m=l  spacecharge  waves 
also  are  Involved,  and  the  complete  quartic  equation  must  be  solved  to  give 
(11). 

For  the  parameters  of  Table  1,  the  low  current  negative  mass  instabil¬ 
ity  has  a  cutoff  at  y  *  7.0.  Since  the  equilibrium  fails  below  y  ~  10.5, 
we  see  that  the  low  current  instability  does  not  arise.  Figure  4,  a  numer¬ 
ical  solution  of  (1)  for  350  A,  illustrates  growth  in  the  low  current 
limit.  The  analytical  result  from  (19)  at  y  *  10.5  exceeds  the  numerical 
by  about  301. 

Rewriting  inequality  (18)  as 


simplifies  comparison  with  (12).  Note  that  (20)  predicts  a  narrow  energy 
bandwidth  at  high  current,  while  (12)  gives  the  opposite.  A  window,  there¬ 
fore,  exists  at  moderate  currents, 

v  (1  +  2zn  a/rb)  -  0. 5-1.0  (21) 

for  which  the  energy  range  of  the  negative  mass  instability  is  minimal. 
For  Instance,  a  750  A,  4.5  MeV  electron  beam  Injected  into  a  betatron  or 
racetrack  accelerator  with  dimensions  as  in  Table  1  exhibits  no  negative 
mass  Instability  whatsoever  for  i  <  6.  Higher  toroidal  modes  are  likely  to 
be  stabilized  by  energy  spread2*4  or  nonlinear  effects.14 


TABLE  1.  NOMINAL  RACETRACK  INDUCTION  ACCELERATOR  PARAMETERS 


Path  Lengths 
Drlfttube  Radius 
Beam  Radius 
Guide  Field 
Beam  Current 
Beam  Energy 
Number  of  Revolutions 
Number  of  Gaps 
Acceleration  per  Gap 
Gap  Resonant  Frequency 
Mode  Quality  Factor 
Gap  Transverse  Impedance 
Gap  Width 


L  ■  460  cm 
R  *  7  cm 
a  *  1  cm 

Bg  *  2  kg  (wc  *  1.173  cm-1) 

I  -  1  kA  (v  -  0.0588) 

U  *  0.4  -  40  MeV  (y  *  1.5  -  80) 

50 

N  -  4 

AU  *  0.2  MeV  (Ay  *  0.4) 

880  MHz 
Q  =  60,  6 
15  ohms 
t  «  5  cm 


(Zl/Q  *  0.5) 
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